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Summary 



This thesis is based on observations performed at the Vacuum Tower Telescope at the 
Observatorio del Teide, Tenerife, Canary Islands. We have used an infrared spectropo- 
larimeter (Tenerife Infrared Polarimeter - TIP) and a Fabry-Perot spectrometer ("Got- 
tingen" Fabry-Perot Interferometer - G-FPI). Observations were obtained during several 
campaigns from 2004 to 2006. We have applied methods to reduce the atmospheric dis- 
tortions both during the observations and afterwards in the case of the G-FPI data using 
image processing techniques. 

We have studied chromospheric dynamics inside the solar disc. The G-FPI provides 
means to obtain very high spatial, spectral and temporal resolution. We observe at several 
wavelengths across the Hor line. With different post-processing techniques, we achieve 
spatial resolutions better than O'.'S. We present results from the comparison of the dif- 
ferent image reconstruction methods. A time series of 55 min duration was taken from 
AR 10875 at ?^ « 36°. From the wealth of structures we selected areas of interest to further 
study in detail some ongoing processes. We apply non-LTE inversion techniques to infer 
physical properties of a recurrent surge. We have studied the occurrence of simultaneous 
sympathetic mini-flares. Using temporal frequency filtering on the time series we observe 
waves along fibrils. We study the implications of their interpretations as wave solutions 
from a hnear approximation of magneto-hydrodynamics. We conlude that a linear theory 
of wave propagation in straight magnetic flux tubes is not sufficient. 

Furthermore, emission above the solar limb is investigated. Using infrared spectro- 
scopic measurements in the Hei 10830 A multiplet we have studied the spicules outside 
solar disc. The analysis shows the variation of the off-limb emission profiles as a func- 
tion of the distance to the visible solar limb. The ratio between the intensities of the blue 
and the red components of this triplet ("R = /biue/^ied) is an observational signature of 
the optical thickness along the light path, which is related to the intensity of the coronal 
irradiation. The observable as a function of the distance to the visible limb is given. 



We have compared the observational with the intensity ratio obtained from Centeno 



( 2006[ ), using detailed radiative transfer calculations in semi-empirical models of the so- 
lar atmosphere assuming spherical geometry . The agreement is purely qualitative. We 
argue that this is a consequence of the limited extension of current models. With the ob- 
servational results as constraints, future models should be extended outwards to reproduce 
our observations. To complete our analysis of spicules we report observational properties 
from high-resolution filtergrams in the Ha spectral line taken with the G-FPI. We find that 
spicules can reach heights of 8 Mm above the limb. We show that spicules outside the 
limb continue as dark fibrils inside the disc. 

One and a half centuries after the hand-drawings by Secchi, the chromosphere is still 
a source of unforeseen and exciting new discoveries. 
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1 Introduction 



This thesis deals with the chromosphere of the Sun. To give some insight to the readers 
which are not familiar with the topics of this work we introduce in Section [TTT] the main 
characteristics of the Sun with a short general description. This will elucidate the position 
of the chromosphere in the solar structure and its role for the outer solar atmosphere. In 
the subsequent Section |1.2[ those aspects of the chromosphere which are treated in the 
present work are specified. Finally Section [L3] indicates the structure of this thesis work. 



1.1 The Sun 



It is just a ball of burning gas 
. . . right? 



The Sun is the central object of the Solar System, which also contains planets and 
many other bodies such as planetoids (small planets), comets, meteoroids and dust parti- 
cles. However, the Sun on its own harbors 99.8% of the total mass of the system, so all 
other objects orbit around it. 

The Sun itself orbits the center of our Galaxy, the Milky Way, with a speed of 217 km/s. 
The period of revolution is ~ 230 million years (the last time the Sun was on this part of 
the Galaxy was the time the Dinosaurs appeared). Compared to the population of stars in 
our galaxy, the Sun is a middle-aged, middle-sized, common type star. In astrophysicist's 
language it is of spectral type G2 and of luminosity class V, located on the main sequence 
of stars in the Hertzsprung-Russell diagram. According to our understandings derived 
from models, it has been on the main sequence for 5 000 million years and it will remain 
there for another 5 000 million years before starting the giant phase. 

The Sun is the closest star to us, the next one being 250 000 times further away, but still 
light from the Sun's surface takes around 8 minutes to reach the Earth. It is the only star 
from where we get enough energy to study its spectrum in great detail and with short tem- 
poral cadence. With indirect methods, we can produce images of the surface structuring 
on other nearby starts. But on the Sun, with current telescopes and techniques, we resolve 
structures down to 100 km size on its surface, which represents approximately the resolu- 
tion limit in this thesis work. We can also investigate the structure of its atmosphere and 
the effects of its magnetism. Actually, we are embedded in the solar wind that has its ori- 
gin in the outer solar atmosphere, the corona of the Sun. Thus, we can make in-situ mea- 
surements. With special techniques and models, we can reconstruct the properties of its 
interior. 
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1 Introduction 




The Sun is the most brilliant object in 
the sky, 12 orders of magnitude brighter 
than the second brightest object, the full 
Moon, which actually only reflects the 
sunlight. Its light warms the surface of the 
Earth and is used by plants to grow. Its ra- 
diation is the input for the climate. The so- 
lar wind separates us from the interstellar 
medium. The magnetism of the Sun pro- 
tects us from cosmic high-energy radiation 
and it influences the climate on Earth. Vio- 
lent events in the solar ultraviolet radiation 
and the solar wind can also disrupt radio 
communications . 

The Sun possesses a complex structure. Essentially, it can be described as a giant 
conglomerate of Hydrogen and Helium (~ 74% and ~ 24% of the mass, respectively) and 
traces of many other chemical elements. Due to its big mass the self-gravitation keeps 
the structure as a sphere. From the weight of the outer spherical gas shells the pressure 
increases towards the center of the sphere. During the gravitational contraction of the 
pre-solar nebula towards its center, i.e. when forming a protostar, the gas has heated up 
by converting potential energy into thermal (kinetic) energy. This produces, together with 
a high gas density, a high pressure, which prevents the sphere from collapsing further 
inwards. Eventually, near the center, temperature and pressure are high enough to ignite 
nuclear reactions. 



Figure 1.1: The apparent size of the Sun on 
the sky is ~ 32', a little bit larger than one 
half degree. 



Structure 

At the core of the Sun the density and temperature (of the order of 13 million Kelvin, 
or 13 • 10^ K) are high enough to fuse hydrogen and burn it into helium. This process 
also produces energy in the form of high-energy photons. This continuous, long-lasting 
energy output from the nuclear reactions keeps the core of the Sun at high temperature 
to sustain the gravitational load from the outer gas shells. Due to the high density, the 
photons are continuously absorbed and re-emitted by nearby ions, and in this way the 
big energy output is slowly radiated outwards, while, towards the surface of the Sun, 
the density decreases exponentially, along with the temperature. Photons reaching today 
the Earth's surface were typically generated on the early times of Homo Sapiens, as the 
typical travel time is ~ 170000 years ( |Mitalas and SiUs|1992[ ). 



At a distance from the center of approximately 70% of the solar radius, the radiation 
process is not efficient enough to transport the huge amount of energy produced in the 
core. There the gas is heated up, and expands, it becomes buoyant and rises. This creates 
convection cells in which hot material is driven up by buoyancy while cool gas sinks to 
the bottom of the cells, where it is heated again. These gas flows transport the energy 
to the outer part of the Sun, where the temperature is measured to be ~ 5 700 K and the 
density is low enough that the photons can escape without much further absorption. The 
outer region from where we receive most of the optical photons can be called the surface 
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NOM 0875, Apr, 26,2006 




Figure 1.2: High resolution image of the "surface" (photosphere) of the Sun with a reso- 
lution of ~ 140 km. Granules are seen all around the photosphere outside the dark areas. 
They form the uppermost layers of the convection zone, in which the energy is transported 
from deep down outwards via gas motions. At the top, the gas cools down by radiating 
photons into space. Localized strong magnetic fields can also emerge and are seen as dark 
areas, the sunspots, which are a consequence of the less efficient energy transport. 



of the Sun, although it is not a layer in the solid state. It is called the photosphere (sphere 
of light). Most of the photons we receive come from this layer are in the visible part of the 
spectrum: light. This is why Nature favored in the late evolution process the development 
of vision instruments that are more sensitive in the spectral region in which most emission 
from the Sun occurs. 

Further out of this layer the atmosphere of the Sun extends radially, with decreasing 
density. In this outer part, with its low density, magnetic fields rooted inside the Sun cease 
to be pushed around by gas flows. This transition occurs together with a still not com- 
pletely understood increase of temperature up to several million degrees. Therefore, there 
must be a layer with a minimum temperature. Standard average models place it at a height 
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of about 500 km with a temperature of about 4000 K, which is low enough to allow the 
formation of molecules like CO or water vapor. Beyond this layer the temperature rises 
outward. Again in standard models, the layer following the temperature minimum has 
an extent of about 1 500 km and its temperature rises to 8 000 - 10000 K. This layer is 
called the chromosphere. The present work deals with some of its properties. Outside 
the chromosphere, the temperature rises abruptly within the transition region. The outer- 
most part of the atmosphere, called corona, drives a permanent outwards flow of particles 
moving along the magnetic field lines. This solar wind extends to 100 000 times the solar 
radius, far beyond Pluto's orbit, to the outer border of our Solar System, the heliopause. 
There the interaction with the interstellar medium creates a shock front, which is being 
measured these years by the Voyager 1 and Voyager 2 probes. 

Beyond this layered structure, the Sun is far more complex. Some other properties, 
which we describe shortly, are: 

- The Sun vibrates. As a self gravitating compressible sphere, it vibrates. Pressure 
and density fluctuations mainly generated by the turbulent convection, are propagated 
through the Sun. Waves with frequencies and wavelengths close to those of the many 
normal modes of vibration of the Sun add up to a characteristic pattern of constructive 
interference. This vibration, although of low amplitude with few 100 m/s in the photo- 
sphere, can be measured and decomposed into eigenmodes by means of Doppler shifts 
and observations of long duration. The propagation of the waves depends on the prop- 
erties of the medium. It is possible then to infer these properties from the measured 
vibration patterns. Some waves propagate only close to the surface, but others can prop- 
agate through the entire Sun. These latter waves provide means to infer some structural 
properties, such as temperature, of the solar interior and test models of the Sun. Global 
Helioseismology provides means to infer the global properties of the interior of the Sun 
studying the vibration pattern, while local helioseismology can depict the surroundings of 
the local perturbations. 

- The Sun rotates. The conservation of angular momentum of a slowly rotating cloud 
that will form a star result, upon contraction, a rapid rotation. It is commonly accepted 
that most of the Sun's angular momentum was removed during the first phases of the 
life of the Sun by braking via magnetic fields anchored in the surrounding interstellar 
medium and by a strong wind. The remaining angular momentum leads to today's solar 
rotation period. But being the Sun not a rigid body this rotation varies from layer to 
layer and with latitude. Gas at the equator rotates at the surface with a period of 27 
days, faster than at the poles where the rotation period is approximately 32 days. Using 
helioseismology observations we know that this diff'erential rotation continues inside the 
Sun, until a certain depth, from which on the inner part rotates like a rigid sphere with 
a period of that at middle latitudes on the surface. This region corresponds to the layer 
where the convection starts, at around 0.7 solar radii, and is called the tachocline. The 
differential rotation creates meridional flows of gas directed towards the poles near the 
surface and towards the equator near the bottom of the convection zone. 

- The Sun shows (complex) magnetic activity. The Sun possesses a very weak overall 
magnetic dipole field. However, the solar surface can host very strong and tremendously 
compUcated magnetic structures, which can be seen through their effects on the solar 
plasma, e.g. less efficient energy transport (that leads to dark sunspots). All matter in 
the Sun is in the form of plasma, due to the high temperature. The high mobility of 
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charges that characterizes the plasma state, makes it highly conductive, causing magnetic 
field lines to be "frozen" into it. Provided that the gas pressure is much higher than the 
magnetic pressure, the magnetic field lines follow generally the dynamics of the plasma. 
The source of these localized strong magnetic fields is still to be understood. The dynamo 
theory addresses this problem suggesting that the weak dipolar magnetic field is amplified 
at the bottom of the convection zone by the stochastic mass motion and shear produced 
by the convection and the difi'erential rotation. 

- The Sun has cycles. The Sun suffers fluctuations in time. Changes occur in the total 
irradiance, in solar wind and in magnetic fields. They happen in approximately regular 
cycles, like the 1 1 years sunspot cycle, and aperiodically over extended times, like the 
Maunder Minimum (a period of 75 years in the XVII century when sunspots were rare, 
and which coincided with the coldest part of the Little Ice Age). These fluctuations mod- 
ulate the structure of the Sun's atmosphere, corona and solar wind, the total irradiance, 
occurrence of flares and coronal mass ejections and also indirectly the flux of incoming 
high-energy cosmic rays. None of these variations are fully understood and their effect on 
the Sun itself or Earth is still under debate. The generally accepted idea about the cyclic 
and more aperiodic fluctuations is that they are caused by variable magnetic fields. These 
are generated by dynamo mechanisms. 

- The Sun evolves. The Sun is now in its main-sequence phase, where the main source 
of energy is the nuclear fusion of hydrogen to helium. After the initial phase of accretion 
of mass, a self gravitating star enters this phase, which lasts for most of its life. In the 
case of the Sun this phase will continue for approximately another five million years, after 
which the later evolution stages include a complex variation of the radius, with burning 
of helium as the source of energy in a later red giant phase. After this stage, the mass of 
the Sun is believed to be not large enough to undergo further fusion stages, and the Sun 
will slowly faint as a white dwarf star. 

Readers can find further general information about the Sun in e.g. |Wikipedial(|Sun|), 



Stix (2002) and many references therein. 



1.2 The chromosphere 

In our short description of the Sun's structure we stated that the atmosphere of the Sun 
comprises a layer above the photosphere in which the temperature begins to rise again 
until the transition region where an abrupt increase of temperature, from approximately 
10000 K to 1 million K, occurs. This first layer above the photosphere is called chromo- 
sphere. The name comes from the greek of "color sphere", as it can be seen as a ring of 
vivid red color around the Sun during total solar eclipse^ 

The boundaries of the chromospheric layer are very rugged, resembling more cloud 
structures than a spheric surface. Above quiet Sun regions the chromosphere can be about 
2 000 km thick, but some structures seen in typical chromospheric lines can reach to much 
higher altitudes, like filaments (that can reach heights of 350000 km). 

The solar chromosphere is a highly dynamic atmospheric layer. At most wavelengths 
in the optical range, it is transparent due to the fact that its density is low, much lower 

'The apparent size of the Sun on the sky happens to be very similar to the apparent size of the Moon, 
leading to annular or total solar eclipses, during which the red ring can be seen. 
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NOAA 0875. Apr, 26,2006 
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Figure 1.3: High resolution filtergram taken in the center of the Ha spectral line, showing 
the chromosphere of the Sun with an image resolution of ~ 150 km. The same field of 



view as image 1.2 The localized strong magnetic fields causing sunspots in the photo- 
sphere are seen now as fibrils around the sunspots. Given the low /3 parameter, the plasma 
is forced to follow the magnetic lines, providing visible tracers and the variety of struc- 
tures seen in the chromosphere. In the image we can see a carpet of spicules, plage region 
and a top view of a rising twisted magnetic flux tube above the active region. This image 
corresponds to the dataset "sigmoid" studied in Chapter |4| 



than in the photosphere below it. Nevertheless, in strong lines like Ha (at 6563 A) or 
Call K and H (at 3934 A and 3969 A, respectively) we have strong absorption (and re- 
emission) which allows direct studies about its peculiar characteristic, like bright plages 
around sunspots, dark filaments across the disk, as well as spicules and prominences 
above the limb. Indeed, recent works, e.g. I Tziotziou et al.| (|2003 ), suggest that many 
of these chromospheric features could all have the same physical properties but within 
different scenarios. 

The temporal evolution of the chromospheric structures is complex. The dynamics of 
a magnetised gas depends on the ratio of the gas pressure Pgas to the magnetic pressure 
Pmag, the plasma /3 parameter, /3 = Pgas/Pmag, with P^^g = B^/(Sn) and B the magnetic 



14 



1 .3 Aim and outline of this work 



field strength q] From the low chromosphere into the extended corona, this plasma pa- 
rameter decreases from values /3 > I, where the magnetic lines follow the motion of the 
plasma (as in the photosphere and solar interior) to a low -beta regime, «; 1, where 
the plasma motions are magnetically driven, and the plasma follows the magnetic field 
lines, creating visible tracers of the magnetism. These eff"ects give rise to a new variety 
of energy transport and phenomena, like magnetic reconnection, filaments standing high 
above the chromosphere or erupting prominences. 



1.3 Aim and outline of this work 



Since the discovery of the chromosphere and since the hand-drawings of |Secchi (1877 1 we 



have been able to observe this solar atmospheric layer in much detail. Many theoretical 
models have been proposed to understand its peculiar characteristics. But, only in the 
last recent years we have been able to address the problem with fine spectropolarimetry 
and high spatial resolution. We can study the fine details and resolve small structures, 
following their dynamics in time. Within these recent advances it has been possible both 
to test current theories and to observe new unexpected phenomena. This work thus aims 
at contributing to the understanding of the solar chromosphere. 

This first Chapter provided a broad introduction to the context of this work. We have 
briefly presented some general properties of the Sun and the chromosphere. In the fol- 
lowing pages, throughout Chapter 2, we summarize some theoretical concepts of radiative 
transfer and spectral line formation needed for this work. We also present general char- 
acteristics of the two spectral lines studied: Ha and Hei 10830 A. Chapter 3 presents in 
detail the observations. There we also summarize the characteristics of the used telescope 
and optical instruments, as well as the data reduction and post-processing methods ap- 
plied to achieve spatial resolutions better than O'.'S. Next, in Chapter 4, we discuss results 
from data on the solar disc, dealing with the chromospheric dynamics and fast events ob- 
served in our data. We present the observations of magnetoacustic waves as well as other 
fast events. Chapter 5 is devoted to the spicules above the solar limb. The analysis of the 
spectroscopic intensity profiles from spicules in the infrared spectral range can be used to 
compare current theoretical models with observations. Further, we present high resolution 
images in Ha of spicules. Finally, the concluding Chapter 6 of this thesis summarizes the 
main conclusions and gives an outlook for future work. 



■^It is very common in astrophysics, specially in solar physics, to use magnetic field strength synony- 
mously with magnetic flux density. The reason is that in most astrophysical plasmas B=H in Gaussian units. 
We follow this use in this thesis. 
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2 Spectral lines 



Most of the information from the extraterrestrial cosmos, also from the Sun, arrives as 
radiation from the sky. It comes encoded in the dependence of the intensity on direction, 
time and wavelength. Also, the polarization state of the light contains information. These 
characteristics of the light we observe from any object have their origin in the interaction 
of atoms and photons under the local properties (temperature, density, magnetic field, 
radiation field itself, . . . ). 

To extract this encoded information from the recorded intensities it is important to un- 
derstand how the radiation is created and transported in the cosmic plasmas and released 
into the almost empty space. 

This Chapter describes in the following sections the basis of radiative transfer and 
spectral line formation. We continue discussing the special properties of the spectral lines 

o 

used in this work: the hydrogen Balmer-or Une (named Ha for short) at 6563 A, and the 
He 1 10830 A multiplet. 



2.1 Radiative transfer and spectral line formation 

Light, consisting of photons, interacts with the gas (of the solar atmosphere, in our case) 
via absorption and emission. Let /^(r, t, Q) be the specific intensity (irradiance) at the 
point r in the atmosphere, at time t, and into direction Q., with = 1. We further denote 
by Kx and as the absorption and emission coefficients, respectively. 
Along a distance ds in the direction Q., the change of is given by 

dh = -Kihds + e^ds , (2. 1) 



or 

^ = -K,h + e,. (2.2) 
ds 

We define also the optical thickness between some points 1 and 2 in the atmosphere 

by 



-r 



dTx = -Kxds ; T3i-T32 = - K^ds, (2.3) 



and the source function of the radiation field as 



S, = -. (2.4) 
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In the solar atmosphere, absorption and emission are usually effected by transitions 
between atomic or molecular energy levels, i.e. by bound-bound, bound-free and free- 
free transitions. If collisions among atoms and with electrons occur much more often 
than the radiative processes, the atmospheric gas attains statistical thermal properties such 
as Maxwellian velocity distributions and the population and ionization ratios according 
to the Boltzamnn and Saha formulae. These properties define locally a temperature T . It 



can be shown (e.g. Chandrasekhar 1960) that in these cases, called Local Thermodynamic 
Equilibrium (LTE), the source function is given by the Planck function or black body 
radiation 

S A varies much more slowly with wavelength than the absorption/emission coefficients 
across a spectral line. Thus, within a spectral line, S a can be considered independent of 
A. 

Generally, LTE does not hold, especially in regions with low densities (thus with 
only few collisions relative to radiation processes) and near the outer boundary of the 
atmosphere from where the radiation can escape into space. The solar chromosphere is a 
typical atmospheric layer where non-LTE applies. In this case, the population densities 
of the atomic levels for a specific transition depend on the detailed processes and routes 
leading to the involved levels. 

Equation 2.2 has the following formal solution 

/(r2) = /(Ti)e-^^'-^2^ + r S{T')e-^^''^'^dT' , (2.6) 

or, for the case when ti ^ oo (optically very thick atmosphere) and T2 = (I(t2 = 0) => 
emergent intensity), then 

Jr*co 
SAr\)e-^^dT\. (2.7) 


A second order expansion of S(t^) leads to the Eddington-Barbier relation 

/,(T, = 0)^5,(T, = 1). (2.8) 

This says that the observed intensity 1^ at a wavelength A is approximately given by the 
source function at optical depth = 1 at this same wavelength. In LTE, the intensity then 
follows the Planck function Ba{T{ta = 1)). 

In spectral lines, the opacity is much increased over the continuum opacity. Since 
the temperature decreases with height in the solar photosphere the intensity in spectral 
lines is decreased, and we probe higher and cooler layers. This explains the formation of 
absorption lines in LTE. 

In non-LTE, when coUisional transitions between atomic levels occur seldom and near 
the outer atmospheric border, photons can escape and are thus lost for the build-up of a 
radiation field in the specific transition. Then, the upper level of the transition becomes 
underpopulated and the source function has decreased below the Planck function at the 
local temperature. It follows that, even for constant temperature atmospheres, a strong 
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absorption line can be observed. 

Outside the solar limb, in the visible spectral range, one observes spectral lines (and 
very weak continua) in emission. In spectral lines, high chromospheric structures are seen 
in front of a dark background. 

2.2 Hydrogen Balmer-of line (Hof) 

Ha at 6563 A is a strong absorption line in the solar spectrum for two reasons: 1) hydro- 
gen is the most abundant element in the Sun, and in the Universe. 2) The Sun, as a G2 V 
star, has the appropriate effective temperature T^ff ^ 5 800 K to have the second level of 
hydrogen populated and thus to make absorptions in Ha possible. 

As all strong lines. Ha possesses a so-called Doppler core and damping wings. The 
Doppler core of Ha and of other Balmer lines is much broader than of other strong lines 
from metals (atomic species with Z > 2). The reason is the large thermal velocity of 
hydrogen compared to that of metals, thus leading to large Doppler widths 

where Aq is the rest wavelength, c the speed of light, 'R the universal gas constant, T 
the temperature and fi the atomic weight (H has the minimum value among the chemical 
elements of fi = 1.008). An eventual "microturbulent" broadening has been omitted in 
Eq.|Z9j 

Another property of the Balmer transitions between the according hydrogen levels is 
the following: Chromospheric lines such as the Can H and K and the Mgi h and k lines 
are weakly coupled to the local temperature through coUisional transitions, effected by 
electrons, between the involved energy levels. Thus, these lines still contain information 
about the temperature of the electrons, although only in a "hidden" manner. 

However, for the Balmer lines of hydrogen and here especially for Ha, there exist 
also the routes for level populations through radiative ionization to the continuum and 
radiative recombination. These routes are taken much more often than the coUisional 
transitions between the involved levels. The ionizing radiation fields, i.e. the Balmer and 
Paschen continua, originate in the lower to middle photosphere and are fairly constant, 
irrespective of the chromospheric dynamics. Only when many high-energy electrons, as 
during a flare, are injected into the chromosphere the Ha line reacts to temperature and 
gets eventually into emission. 

Nonetheless, the chromosphere observed in Ha exhibits rich structuring, due to ab- 
sorption by gas ejecta, due to Doppler shifts of the Ha profile in fast gas flows along 
magnetic fields, and due to channeling of photons around absorbing features. 

2.3 He 1 10830 A multiplet 

Helium is the second most abundant element in the Universe, also in the Sun. It was 
first discovered in the Sun in 1868 (from where it was named after the greek word of 
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Sun). At the typical chromospheric temperatures there is not enough energy to excite 
electrons to populate the upper levels from where these transitions occur. In coronal holes 
the helium lines are substantially weaker compared with the quiet Sun outside the limb. 
More information about recent advanc es in measur ing chromospheric magnetic fields in 
the He I 10830 A line can be found in 



Lagg(2007). 



The energy levels that take part in the transitions of the Hei 10830 A multiplet are 
basically populated via an ionization-recombination process ( [Avrett et al.||1994| ). The 
much hotter corona irradiates at high energies both outwards to space and inwards, to 
the chromosphere. The EUV coronal irradiation (CI) at wavelengths lambda A < 504 A 
ionizes the neutral helium, and subsequent recombinations of singly ionized helium with 
free electrons lead to an overpopulation of the upper levels of the Hei 10830 multiplet. 

Alternative theories suggest other mechanisms that may also contribute to the forma- 
tion of the helium lines via the colhsional excitation of the electrons in regions with higher 
temperature (e.g. [Andretta and Jones|1997[ ). 
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The Hei 10830 A multiplet consists 
of the three transitions of the orthohelium 
(total spin of the electrons S = l) energy 
levels, from the upper term with angu- 
lar momentum L = 1 to the lower with 
L = 0, in particular from ^P2,i,o, which has 
three sublevels (7 = 2, 1,0), to the lower 
metastable term ^Si, which has one single 
level (7=1) (see Fig. |2.1[ ). The two tran- 
sitions from the J=2 and J=l upper levels 
appear mutually blended, i.e. as merely 
one line, at typical chromospheric temper- 
atures, and form the so-called red compo- 
nent, at 10830.3 A. The two red transi- 
tions are only 0.09 A apart. The blue com- 
ponent, at 10829.1 A, corresponds to the 
transition from the upper level with J=0 to 
the lower level with J = 1 . 

The formation height of these lines is believed to be between 1 500 and 2 000 km, (e.g. 
Centeno|2006 1 although, as we already mentioned, the chromosphere is strongly rugged. 



red component 
10830.3 

J=l 



Figure 2.1: Schematic Grotrian diagram for 
the Hei 10830 A multiplet emission lines. 



The Lande factors of the lines are not zero, meaning that they are sensitive to external 
magnetic fields. 

A more detailed description about the properties of the Hei 10830 multiplet, in par- 
ticular related to the emission profiles observed in spicules above the limb is given in 
Chapter [5} 
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For the present work we used data from two different instruments, both mounted on the 
same telescope, the Vacuum Tower Telescope (VTT, Sec. 3.2) in Tenerife. One of the in- 
struments, the Gottingen Fabry-Perot Interferometer (G-FPI, Sec. 3.3.1 ) is able to achieve 
very high spatial resolution while the other, the Tenerife Infrared Polarimeter (TIP, Sec. 



3.4.2), is able to obtain full Stokes spectropolarimetric data with very high spectral res- 



olution. Both instruments, in combination with the Kiepenheuer Adaptive Optics System 
(KAOS, Sec. |3.2.1[ ), provided the data for this work. 

In this Chapter we will describe the telescope, the instrumentation, the observations, 
and the reduction techniques. The latter are aimed at removing as many instrumental 
effects as possible. 



3.1 Angular resolution and Seeing 

When using any kind of an optical imaging system, the angular resolution in the focal 
plane is limited by diffraction at the aperture of the instrument. For circular apertures the 
image of a point source (the PSF) is an Airy function with a certain Full Width at Half 
Maximum (FWHM). Two point sources closer than the FWHM of a certain instrumental 
PSF are difficult to distinguish. If one considers diffraction of a telescope with a circular 
aperture of diameter D, the angular resolution limit is, in the usual Rayleigh definition, 

Qmin = 1-22 ^ . (3.1) 

The factor 1 .22 is approximately the first zero divided by n of the Bessel function involved 
in the Airy function. In the focal plane of such a telescope with a focal length / the spatial 
resolution is therefore d = 1,22 A f/D. For good sampling this should correspond to, or 
even be larger than, the resolution element of the detector (2 pixels). In the case of the 
VTT, with a main mirror of D = 70 cm, the diffraction limited resolution is Ql.'24 at 6563 
A (Ha) and (y.'39 at 10830 A (Hei triplet). In solar observation, it is common to use as 
the diffraction limit simply a^in = A/D. At this angular distance the modulation transfer 
function (MTF) has become zero. 

Unfortunately all imaging systems on the ground are subject to aberrations that de- 
grade the image quality, resulting in a much lower spatial resolution than the diffraction 
limit. The light we observe from the Sun travels unperturbed along approximately 150 
million km, but during the last few microseconds before detection it becomes distorted 
due to its interaction with the Earth's atmosphere and our optical instrument. 

The refraction index of the air is very close to 1 at optical wavelengths, but depends 
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on the local pressure and temperature. Their fluctuations in space and time produce aber- 
rations of the wavefronts from the object to be observe(Q Since the time scale of the 
variation of the aberrations ofx 10 ms is usually smaller than the integration time, it also 
produces smoothing of the image details. Thus, the information at small scales is lost. 

Further, the turbulent state of the air masses through which the light is passing varies 
on small angular scales. This produces an anisoplanatism of the wavefronts arriving at 
the telescope, with angular sizes of the isoplanatic patches not larger than approximately 
10". 

Beside the atmospheric factors, the final quality of the image is influenced by local 
factors like the aerodynamical shape of the telescope building or convection around the 
building and the dome. 

Finally, the internal seeing of the telescope plays an important role for the image qual- 
ity. Convection along the light path in the telescope triggered by heated optical surfaces 
can be avoided by allowing air flowing freely through the structure or, quite the contrary, 
by evacuating the telescope. 

In solar physics we usually measure the average image quality of the observations 
estimating the diameter of a telescope that would produce, from a point source, an image 
with the same diflraction-limited FWHM as the atmospheric turbulence or internal seeing 
would allow even with a much larger telescope aperture. This is called the Fried parameter 
(ro). Typically, upper limits for the "Observatorio del Teide" are ro « 15 cm during night- 
time and ro « 7 cm during day. 

Besides these structural requirements for best seeing conditions, there are nowadays 
methods for correcting the images for seeing distortions to obtain near diffraction limited 
resolution. In this thesis we have used various methods: We correct partially the aberra- 
tions in real time using adaptive optics (Sec. |3.2.1[ ) which can increase the ro around the 
center of the field of view up to ro ~ 25 cm and we also apply post-processing methods of 
image reconstruction (Sec. |3.3.3[) to approach the upper limit of ro < 70 cm. 



The local values of the temperature and pressure depend on the complicated turbulent dynamics of 
the atmosphere. This includes friction and heating of the Earth's irregular surfaces, condensations and 
formation of clouds, shears produced by strong winds, . . . For more information we refer to e.g. Saha 2002 
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3.2 Telescope 



The Vacuum Tower Telescope (VTT, |Soltau|[T985t Fig. |3.1| ) is located at the Spanish 
"Observatorio del Teide" (2400 m above sea level, 16°30' W, 28° 18' N) in Tenerife, Ca- 
nary Islands. It is operated by the Kiepenheuer-Institut fiir Sonnenphysik, Freiburg, with 
contributions from the Institut fiir Astrophysik in Gottingen, the Max-Planck-Institut for 
Sonnensystemforschung, Katlenburg-Lindau, and the Astrophysikalisches Institut Pots- 
dam. 



The VTT optical setup is depicted in Fig. [X2 
At the top platform of the building, a coelostat 
achieves to follow the path of the Sun on the sky, 
by means of two flat mirrors of very high optical 
quality. The primary coelostat mirror rotates clock- 
wise (seen pole-on) about an axis which is con- 
tained in the mirror surface and is parallel to the 
Earth's rotation axis. It reflects the sunlight towards 
the secondary mirror. The latter redirects the beam 
towards the fixed telescope in the tower. The tele- 
scope is an ofl'-axis system. It consist of a slightly 
aspherical main mirror of 70 cm diameter and a 
focal length of 46 m, and of a folding flat mirror. 
The free aperture of the circular entrance pupil with 
D=70 cm gives the telescopic diffraction limit for 
the angular resolution of a„„„ = AjD « 0'.'16 for A 
in the visible spectral range. 




Figure 3.1: Building which houses 
the solar Vacuum Tower Telescope. 



To avoid turbulent air flows inside the telescope 
caused by heated surfaces, the telescope is mounted 

in a tank that is evacuated to 1 mbar. The vacuum tank has high quality transparent 
entrance and exit windows located below the coelostat and close to the primary focus, 
respectively. 

Shortly after the entrance window, a small part of the sunlight is reflected out to a 
second imaging device. This uses a quadrant cell to track the image of the solar disc and 
to correct slow image motions, e.g. due to a non-perfect hour drive of the coelostat. Tele- 
scope pointing to a target inside and near the solar disc is achieved by moving this tracking 
device as a whole in the image plane. The imbalanced illumination of the quadrant cell is 
transformed to a tip-tilt motion of the secondary coelostat mirror. 



After the main vacuum tank, the adaptive optics (Sec 3.2.1 1 device is located. This 
optical system is able to correct in real time the low order aberrations of the incoming 
wavefronts of the light beam caused by the turbulence in the Earth's atmosphere. After 
the adaptive optics system, which can optionally be moved in or out of the path, the light 
path continues to the vertical slit spectrograph or to a folding mirror that can be used to 
direct the light to diff"erent other available science instruments. 
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Figure 3.2: Optical setup of the VTT. The coelostat (mirrors ml,m2) follows the path of 
the Sun on the sky and directs the light to the entrance window of the vacuum tank (blue 
shaded). Mirror m3 takes out a small amount of the light and feeds the guiding telescope 
mounted outside the vacuum tank. The coUimating mirror m5 produces, together with the 
flat mirror m(5, the solar image in the primary focal plane behind the exit window of the 
vacuum tank. There, a flat mirror can be mounted under 45° to the vertical (not shown) to 
feed post-focus instruments in optical laboratories. The adaptive optics system is located 
below the exit window, and it is used optionally. 



3.2.1 Kiepenheuer Adaptive Optics System 



As mentioned in the beginning of this Chapter (Section 3.1 1 the atmosphere of the Earth 
degrades the quality of the images during observations. KAOS (Kiepenheuer Adaptive 
Optics System, von der Luhe et al.|2003[|Berkefeld|2007 1 is a realtime correction device 



that calculates and corrects the instantaneous aberrations of the wavefront using special 
deformable mirrors. 



The optical scheme of a typical adaptive optics (AO) system is shown in Fig. 3.3 By 
means of a dichroic semitransparent beam splitter, part of the light entering the system 
is directed to the wavefront sensor. The latter, a Shack-Hartmann sensor, consists of a 
lenslet array positioned in an image of the entrance pupil and a fast CCD detector. Each 
lenslet, cutting out a subaperture of the pupil image, produces an image of a small area 
on the Sun on a subarea of the CCD. Using a good, i.e. as sharp as possible, subimage of 
the present scenery on the Sun and with a correlation algorithm, it is possible to compute 
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Figure 3.3: Scheme of typical AO. Inside the closed loop, a fraction of the incoming light 
is directed to the KAOS camera (semitransparent mirror ml), where a lenslet array (//) 
produces many subfield images with light from different parts of the pupil. The calculated 
instantaneous aberration is compensated using the two (tip&tilt and deformable) mirrors, 
every 0.4 ms. 



the displacement of each subimage and to estimate from this the aberrations of the wave- 
front. Every aberration can be expressed by a sum of adequate polynomials (for example 
Zernike polynomials) with appropriate coefficients. Each polynomial represents a specific 
wavefront aberration, e.g. tilt, defocus, astigmatismus . . . The AO is able to correct the 
low orders of the aberration, that is those with the largest scales. For this purpose it has 
two active optical surfaces (both of them in the main lightbeam, so the correction is done 
in a closed loop). In the case of KAOS the first element is the tip-tilt mirror that is able to 
displace the whole image in two perpendicular directions, thus tracking on the reference 
image. The second optical element is a bymorphous deformable mirror with 35 actuators. 
With appropriate voltages, the surface of this mirror obtains a shape that corrects the aber- 
rations of the incoming wavefront up to the 27*'' Zernike polynomial. This correction is 
done in a fast closed loop at 2100 Hz. The bandwidth of KAOS is 100 Hz. It thus operates 
at timescales comparable to that of the variation of the turbulence in the atmosphere. 

As already mentioned, the aberration of the wavefront is not constant, i.e. not isopla- 
natic across the whole field of view (FoV). The wavefront camera has a restricted FoV 
of 12" X 12" where the assumption of isoplanatism is approximately valid. The center 
of this subfield of AO correction is called lockpoint. The restricted area of isoplanatism 
is one of the main limitations of current AO systems. The corrections are calculated for 
the lockpoint feature we are tracking on and applied to the whole FoV of the telescope. 
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Therefore the correction becomes increasingly inaccurate with increasing distance from 
the lockpoint. The quality of the image is degraded outwards from the center of the FoV, 
where the lockpoint is usually located. Fortunately this can be taken into account using 
post factum image reconstruction like speckle interferometry and blind deconvolution. 

In night-time astronomy, AO systems lock on a star image, so the displacements of 
the subfields imaged by the lenslet are easily calculated. In solar observations, the im- 
age used by the AO comes always from an extended source, making the calculations of 
the displacements much more demanding. In solar AOs, a reference image is taken and 
updated regularly during operation, and correlations between this image and the subfield 
images are used. For well defined maxima of the correlation functions we need features 
with sufficient contrast inside the FoV to lock on with the algorithm, e.g. a pore or the 
granulation pattern. Moreover, the wavefront sensor can only work with a high light level, 
e.g. integrated over some wavelength. So it is not possible to lock for example on fea- 
tures within the Hot line with low intensity. Also, as we will explain in Sec. |3.4.2[ near or 



off-limb observations are difficult as the AO algorithm is not able to track on that kind of 
references, as the one-dimensional limb image. 



3.3 High spatial resolution 

For our study of the dynamics of chromospheric structures, we are interested in observa- 
tions with the highest possible spatial resolutiorj^ with the highest achievable temporal 
cadence, and with as much spectral information as possible. For that purpose we used the 
"Gottingen" Fabry-Perot Interferometer (G-FPI). Here, the designation FPI stands as pars 
pro toto, for the whole post-focus instrument, a two-dimensional spectrometer based on 
wavelength scanning Fabry-Perot etalons. It was developed at the Universitats-Stemwarte 
Gottingen ( |Bendlin et al.fT992l|Bendlin| 19931 [Bendlin and Volkmer|1995| ). Subsequently, 



it had undergone several upgrades ( [Koschinsky et al.||2001[ [Puschmann et al.|2006[ [Bello 



Gonzalez and Kneer||2008j ). For the present work, the G-FPI with the high-efficiency 



performance described by |Puschmann et aT ( 2006| ) was employed. 



Basically, this instrument was able, at the time the data for this study were taken, to 
produce an image from a selected wavelength range with a narrow passband of 45 mA 
FWHM at 6563 A (Ra). A recent upgrade has reduced the FWHM. The spectrometer 
also can be tuned to almost any desired wavelength, being able to scan a spectral line, 
producing 2D filtergrams (images) at, e.g., 20 spectral position along a line. If we scan 
iteratively one spectral line we obtain a time sequence of very high spatial resolution, at 
several spectral positions and with a cadence which would be the time required to scan 
the full line, which is typically in the order of 20 seconds for our data. 

The main limitation of this kind of observational procedure is that the images corre- 
sponding to a single scan are not obtained simultaneously, as they are taken consecutively. 
This is of special importance when we compare the images in the two wings of a spec- 
tral line, as the small-scale solar structure under study may have changed during the time 
needed to scan between these positions. This should be taken into account when studying 
features whose typical timescale of variation is comparable to the scanning time. In Sec. 



^It has become a widespread custom in solar observations to use "spatial resolution" synonymously with 
"angular resolution". 
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4.2.1[ we will see that this limitation can partly be compensated when we have a long 



temporal series. 
3.3.1 Instrument 



The Gottingen Fabry-Perot Interferometer ( 


Bendlin and Volkmer 1995[ Volkmer et al. 


19951 Koschinsky et al.|2001 


Puschmann et al. 12006]) is a speckle-ready two-dimensional 



(2D) spectrometer. It is able to scan a spectral line producing a set of speckle images 
at several spectral position with a narrow spectral FWHM, while taking simultaneous 
broadband images, needed for the post factum image reconstruction. 

Fabry-Perot interferometer (FPI) 

A Fabry-Perot interferometer, or etalon, is an interference filter possessing two plane- 
parallel high-reflectance layers of high quality (~ /I/ 100). Light entering the filter is 
many times reflected between the plane-parallel reflecting surfaces. These reflections will 
produce destructive interference for transmitted light at all wavelengths but the ones for 
which two times the spacing d of the plates is very close to a multiple of the wavelength. 
This eff"ect gives rise to a final Airy intensity function ( Born and Wolf|1999 ): 



I — Imax 4^ -25' (3-2) 

where the maximum intensity Imax = (ir^ > T is the transmittance, R is the reflectance 
{R = 1 - r if absorption is negligible), and the dependence on wavelength A, angle of 
incidence 0, and refractive index n of the material between the surfaces is 

An 

6 = —nd COS©. (3.3) 
A 

The narrow transmittance of the filter can be tuned to any desired wavelength by 
changing the spacing d (or the refractive index n, for pressure controlled FPIs). One 
single FPI produces a channel spectrum according to the interference condition, i.e. for 
normal incidence (0 = 0°) and assuming n=l, 

mA = 2d (3.4) 

with m being the order. From here, the distance to the next transmission peak, or free 
spectral range (FSR), follows as 

FSR = ^. (3.5) 

To suppress all but the desired transmission, the G-FPI has a second Fabry-Perot etalon 
with different spacing, i.e. different FSR. Both Fabry-Perot etalons need to be synchro- 
nized when scanning in order to keep the desired central transmittance peaks coinciding. 
The combination of two FPI with different FSR removes effectively the undesired trans- 
mission peaks from other orders. An additional interference filter (FWHM 8 A) is used 
to reduce the incoming spectral range to the spectral line under observation. The combi- 
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Figure 3.4: Example of the narrow-band scanning with the G-FPI. Left: One narrow- 
band frame from a two-dimensional spectrometric scan through the hydrogen Balmer-or 
line (Ha). Right: Ha line; solid black from the Fourier Transform Spectrometer (FTS) 



atlas (Brault & Neckel, quoted by Neckel 1999); blue: FTS profile convolved with the 



Airy transmission function of the FPIs; dashed average Ha profile observed with the 
spectrometer at 21 wavelength position (rhombi) with steps of 100 mA. The red line is 
the Airy transmission function, positioned at the wavelength in which the image in the 
left panel was taken, and re-normalized to fit on the plot.. 



nation of these three elements produces a single narrow central peak, as depicted in Fig. 



The FP etalons are mounted close to an image of the telescope's entrance pupil in the 
coUimated, i.e. parallel, beam. On the one hand, this avoids the "orange peel" pattern 
in the images, which one obtains with the telecentric mounting near the focus and which 
arises from tiny imperfections of the etalon surfaces. On the other hand, in the coUimated 
mounting one has to deal with the fact that the wavelength position of the maximum 
transmission depends on the position in the FoV. This can be seen from Eq. |3. 3 [ where the 
angle of incidence changes with position in the FoV. 

For the post factum image reconstruction (Sec. 3.3.3) we have to acquire simulta- 
neously short-exposure images from the narrow-band FPI spectrometer and broadband 
images. The latter are taken through a broadband interference filter (FWHM x 50 A) at 
wavelength close to the one observed with the spectrometer. Two CCD detectors, one for 
each channel, with high sensitivity and high frame rates were used which allow a high 
cadence of short exposures. All processes (simultaneous exposures, synchronous FPI 
scanning and observation parameters) are controled by a central computer. The imaging 
on the two CCDs is aligned with special mountings and adjusted to have the same image 
scale on the two detectors. 



The optical setup is shown schematically in Fig. 3.6 From the focal plane following 
KAOS the image from the region of interest on the Sun is transferred via a 1 : 1 re- 
imaging system into the optical laboratory housing the FPI spectrometer. In front of 
the focus at the spectrometer entrance, a beam splitter directs 5% of the light into the 
broadband channel. The latter contains a focusing lens, the broadband interference filter 
(IFl), a filter blocking the infrared light (KGl, from Kaltglas = "cold glass", notation 
by Schott AG), a neutral density filter to reduce the broadband light level, and a detector 
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Figure 3.5: Transmission functions for the narrow -band channel of the G-FPI with the Ha 
setup. The periodic Airy function of the narrow-band FPI (dashed line) coincides in the 
central wavelength with that of the broadband FPI (strong dashed green line). The global 
transmission of both FPIs has one single strong and narrow peak at the central wavelength 
(purple strong line). An additional interference filter (red line) is mounted to restrict the 
light to the scanned spectral line. 



CCDl. 

Most of the light (95 %), enters the narrow-band channel of the spectrometer through 
a field stop at the entrance focus. After the field stop follow: an infrared blocking filter 
(KG2), the narrow interference filter (IF2), a coUimating lens giving parallel light, the two 
Fabry-Perot etalons (FPI-B and FPI-N), a camera lens focusing the light on the detector 
CCD2. Figure [34] gives an example of the type of observation one can obtain with this 
narrow-band spectrometer. 

The instrument has additional devices for calibration and adjustment: a feed of laser 
light, facilities to measure with a photomultiplier and to aid identifying the spectral line 
to be observed, and a feed of continuum light for various purposes, e.g. co-aligning the 
transmission maxima of the etalons or measuring the transmission curve of the pre-filter 
IF2. 



3.3.2 Observations 

For the study of the chromospheric dynamics on the basis of high resolution observations 
we have used three data sets. Table 13. ll lists the details for each data set: 

• Dataset mosaic focuses on the study of a large active solar region, where we find fast 



moving dark clouds, as we will discuss in Sec. 4. 1 These data were obtained before 
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Figure 3.6: Schema of the "Gottingen" Fabry-Perot interferometer optical setup. After 
KAOS, the light is transferred from the telescope's primary focus to the spectrometer. A 
beam splitter BS directs 5% of the light into the broadband channel consisting of a focus- 
ing lens LI, a broadband interference filter IFl (FWHM ^ 50 A), an infrared blocking fil- 
ter KGl ("Kaltglas"), a neutral density filter ND, and the CCDl detector. 95% of the light 
enter the spectrometer through a field stop at the entrance focus. Then follow: infrared 
blocking filter KG2, interference (pre-) filter IFII (FWHM ^ 6A . . . lOA, depending on 
the spectral line and wavelength range), coUimating lens L2, the two FPI etalons FPI B 
and FPI N (FWHM « 45mA at Ha), the focusing camera lens L3 and the CCD2 detector. 
CCDl and CCD2 take short-exposure (3-20 ms) images strictly simultaneously. 



the instrument upgrading in 2005 (Puschmann et al. 2006) with the old cameras. 
The exposure time was six times longer than with the new CCDs and the FoV of a 
single frame is one fourth of that of the new version of the G-FPI. The observers of 
these data were Monica Sanchez Cuberes, Klaus Puschmann and Franz Kneer. 



• Dataset sigmoid uses the improvements of the instrument from 2005 and was ob- 
tained during excellent seeing conditions from a very active region. During the time 
span of our observations at least one flare was recorded from this region in our FoV. 
Our focus with these data is the study of fast events and magnetoacustic waves (Sec. 



4.2.4) with the original intention to detect Alfven waves. Examples of these data 



were also used to compare the results from difl'erent methods of post factum image 



reconstruction, as we will show in Sec. 4.3 



With dataset limb and in Sec. 4.3 we apply blind deconvolution methods for image 



reconstruction (see Sec. 3.3.3.0. The observations were taken with the G-FPI, 



renewed in 2005, to study with very high spatial resolution the evolution of spicules 
as seen in the Ha line. 
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3.3.3 Data reduction 

After the recording of the data, several processing steps have to be carried out in order to 
minimize the instrumental effects. These are mainly to take into account the differential 
sensitivity of the CCDs from one pixel to another or the fixed imperfections on the optical 
surfaces positioned close to one of the focal planes. This concerns for example dust on 
the beam splitter, on the infrared blocking filters and interference filters and the CCDs. 
In this step we also remove an imposed bias signal applied electronically to every frame. 
This is the usual treatment of any CCD data. 

For this purpose we take fiat fields, dark, continuum and target images (see Fig. |3.7[ ). 

Target. A target grid is located in front of the instrument, in the primary focal plane. 
Target frames therefore display in both channels a grid of lines that are used to focus and 
align the cameras in both channels. This is crucial for the image reconstruction. 

Continuum data are taken with the same scanning parameters as with sunlight but 
using a continuum source, so we can test the transmission of the scanning narrow-band 
channel. 

Dark frames are taken with the same integration time but blocking the incident light. 
These frames have information of the differential and total response of the CCD array 
without light, in order to remove this effect from the scientific data. 

Flat fields are frames with the same scanning parameters and with sunlight, but with- 
out solar structures. In this way we can see the imperfections and dust on the optical 
surfaces fixed on every frame taken with the instrument, and remove them dividing our 
science data by these fiat frames. To avoid signatures from solar structures in the flat 
frames, the telescope pointing is driven to make a random path around the center of the 
solar disc far from active regions. 



Data set name 


"mosaic " 


"sigmoid" 


"limb" 


Date 


May,31-^',2004 


April,26"',2006 


MayA"\ 2005 


Object 


AR0621 


AR10875 


limb 


Heliocentric angle 


H = 0.68 


lu = 0.59 


// = 


Scans # 


5 


157 


5 


Cadence 


45 s 


~ 22 s (see Sec. 4.2.1) 


~ 19 s 


Time span 


4 min 


55 min 


2 min 


Line positions # 


18 


21 


22 


FWHM 


50 A broadband / 45 mA narrow-band 


Broadband filter 


6300 A 


Stepwidth 


125 mA 


100mA 


93 mA 


Exposure time 


30 ms 


5 ms 


Seeing condition 


good 


ro « 32 cm 


ro « 20 cm 


KAOS support 


yes 


Image reconstruction 


speckle 


AO ready speckle 


MFMOBD 


Field of view 


33"x 23"(total 103"x94") 


77"x 58" 



Table 3.1: Characteristics of the data sets taken with the G-FPI used in this work. 
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(a) Broad band raw frame 



(b) Flat field frame 






(c) Dark frame 



(d) Reduced frame 



Figure 3.7: Example of the standard data reduction process. Every frame taken with 
the CCD (a) includes instrumental artifacts like shadows from dust particles on the CCD 
chips or the filters near the focus (Fig. b) and the intrinsic diff'erential response of each 
pixel (c). Subtracting the dark frame and dividing by the flat response provides a clean 
frame (d). 



Thus, to reduce the instrumental effects we use the following formula, for each chan- 
nel and for each spectral position independently: 

raw frame - mean dark 

reduced J rame = — — — — ; — - . (3.6) 

mean flat field - mean dark 

Our instruments produce data sets that can be subject to post factum image recon- 
struction. We have appUed speckle and blind deconvolution methods to minimize the 
wavefront aberrations and to achieve spatial resolution close to the difiraction limit im- 
posed by the aperture of the telescope. 

The aberrations are changing in time and space. In a long exposure image, the tem- 
poral dependence will produce the summation of different aberrations, blurring the small 
details of the image. Therefore, for post-processing, all image reconstruction methods 
need input speckle frames with integration times shorter than the typical timescale of 
the atmospheric turbulence. With this condition fulfilled, the images appear distorted 
and speckled but not blurred, and still contain the information on small-scale structures. 
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Another common characteristic of speckle methods is the way to address the field depen- 
dence of the aberrations. In a wide FoV each part of the frame is affected by different 
turbulences. That is, inside the atmospheric column affecting the image, there are spatial 
changes of the wavefront aberration. Therefore, the FoV is divided into a set of overlap- 
ping subfields smaller than the typical angular scale of change of the aberrations (5"- 8"), 
the isoplanatic patch. 

Speckle interferometry denotes the interference of parts of a wavefront from different 
sub-apertures of a telescope. This results in a speckled image of a point source, e.g. 
of a star. The effect is used for "speckle interferometric" techniques of postproccesing. 
They are able to remove the atmospheric aberrations of the wavefronts that degrade the 
quality of the images. In the following Sections we introduce the basic background of the 
methods used and provide some examples and further reference. 



3.3.3.a Speckle interferometry of the broadband images 



This method is based on a statistical approach to deduce the influence of the atmosphere. 
It was developed following the ideas of |Fnedl ( fT965| ), |Labeyrie| ([1970]), |Korff| ([1973]), 
Weigelt] ( |1977| ), |von der Luhe| ( |1984| ) . The code used for our data was developed at 
the Universitats-Sternwarte Gottingen ( |de Boer||1996| ) . The sigmoid dataset uses the lat- 
est improvements to take into account the field dependence of the correction from the AO 
systems ( Puschmann and Sailer|2006 ). 

In what follows we present a brief overview of the method: The observed image (z) 
is the convolution (★) of the true object (o) with the Point Spread function (PS F). The 
PS F is the intensity distribution in the image plane from a point source with intensity 
normalized to one, i.e. 



PSF(x,y)dxdy = 1 



(3.7) 



where the integration is carried out in the image plane. The PS F depends on space, time 
and wavelength. Its Fourier transform (^) is the OTF, Optical Transfer Function 



^(i) = ^(oirPSF) 



1 = 0- OTF . 



(3.8) 



A normal long exposure image would be just the summation of N speckle images: 

N N 



J]h = O.J]OTF, 



(3.9) 



!=1 



(=1 



The OTFi are continuously changing in time, which leads to a loss of information. The 
temporal phase change of the OTFi will, upon this summation, reduce strongly or even 
cancel the complex amplitudes at high wavenumbers. Labeyrie ( 1970[ ) proposed to use 
the square modulus, to avoid cancellations: 



. N N 



(3.10) 



(=1 (=1 

Yet this procedure also removes the phase information on o. Thus, the phases have to be 
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retrieved afterwards. STF is the Speckle Transfer Function, it contains the information on 
the wavefront aberrations during N speckle images. To deduce this STF is therefore one 
of the aims of the speckle method. On the Sun, point sources do not exist. It is thus not 
a trivial task to determine the S TF. There are, however, models of S TF for extended 
sources from the notion that they depend only on the seeing conditions, through the Fried 
parameter ro ( Korff|1973 1. This parameter can be calculated statistically using the spectral 



ratio method ( |von der Luhe||1984| ). As this is a statistical approach, a minimum number 



of speckle frames must be used, more than 100. 

To recover the phases of the original object the code uses the speckle masking method 
( Weigelt||1977[ Weigelt and Wimitzer|[1983j ). It recursively recovers the phases from low 



to high wavenumbers. 

Finally a noise filter is applied, zeroing all the amplitudes at wavenumbers higher than 
a certain value, which depends on the quality of the data. 

Influence of the AO on the speckle interferometry 

As explained in Sec. |3.2.1| the AO systems provide a realtime correction of the low or- 
der aberrations (up to a certain order of Zernike polynomials). Nonetheless, given the 
anisoplanatism of the large field of view, the corrections are calculated for the lock point 
and applied to the whole frame, resulting in a degradation of the image correction from 
the lock point outwards. The problem arises from the different atmospheric columns tra- 
versed by the light from different parts in the FoV. This creates, after the AO correction, 
an annular dependence of the correction about the lock point and therefore an annular de- 



pendence of the STFs when processing the data. Puschmann and Sailer (2006) provided 



a modified version of the reconstruction code that computes different S TFs for annular 
regions around the lock point, providing a more accurate treatment over the field of view. 

The sigmoid dataset was reduced using this last version of the code, improving sub- 
stantially the quality of the results. Both AO and speckle interferometry work best with 
good seeing, and this data set was recorded under very good seeing conditions. 

3.3.3.b Speckle reconstruction of the narrow-band images 

The narrow-band channel scans the selected spectral line, taking several (~ 20) images 
per spectral position. The statistical approach as for the broadband data can not be ap- 
plied given the low number of frames per spectral position. To reconstruct these images 



from this channel we use a method proposed by Keller and von der Liihe ( 1992 1 and im 



plemented in the code by Janssen (2003 1. For each narrow-band frame, there is a frame 
taken simultaneously in the broadband channel, which is degraded by the same wave aber- 
rations. The images in the broadband channel were taken at 6300 A, i.e. at a wavelength 
260 A shorter than that of Ha. We neglect the wavelength dependence of the aberration. 
For each position in the spectral line, for each subfield, we have a set of pairs of simul- 
taneous speckle images from the narrow- and broadband channel, with a common OTFi 
for each realization in both channels: 

= OBroad-OTFi (3.11) 
ONarro..-OTFi (3.12) 
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(a) Average of 330 speckle images (total exposure 
time ~ 1,6 s). 



(b) Single speckle frame, 5 ms exposure time. 




(c) Reconstructed broadband image, using 330 speckle frames. 

Figure 3.8: Example of improvement of broadband images with the speckle reconstruc- 
tion. The size of the image is ~ 34" x 19". The achieved spatial resolution is close to the 
diffraction hmit, 0'.'22, with the diffraction limit a^m = A/D = (y.'\9 at A = 6563 A (He) 
and telescope aperture D = 10 cm. 



Using Equation 3.11 in 3.12 the reconstructed narrow-band image ONamm is obtained 



from the minimization of the error metric 



^ — / J Narrow ' q 
i=l 



Broadi 



-I 



Broadj 



Narrowj 



(3.13) 



where N is the number of images taken at one wavelength position. Minimization of E 
with respect to O^arrow yields 



O 



Narrow 



H 



,•=1 ^Narrowi ' ^Broadj 



o 



Broadi 



(3.14) 
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spatial frequency (1/Mm) 

0.0 1.3 2.6 3.9 5.3 6.6 




54 1.072 0.541 0.362 0.272 0.218 

Nyquist 

angular scale (arcsec) limit 



Figure 3.9: Power spectra showing the influence of the post factum reconstruction. Ordi- 
nate is the relative power on logarithmic scale, and abscissa is the spatial frequency, from 
the largest scales near the origin to the smallest scales at the Nyquist limit, correspond- 
ing to two pixels. A long exposure image {black dotted line), taking the average of all 
speckle images, has very low noise, but the power is also low at all frequencies ^ 0.8 
Mm"' (blurring effect). A single speckle frame {dashed blue line) has more power at all 
frequencies, but also much more noise (more than two order of magnitude). The speckle 
reconstructed frame {red solid line) keeps the noise low while it possesses higher power 
at all frequencies, where the spatial information on small-scale structures is stored. 



Here we have included a noise noise filter {H) to remove the power at spatial frequencies 
higher than a certain threshold above which the noise dominates. The noise power is 
obtained from the flat field data. 



3.3.3.C Multi object multi frame blind deconvolution (MOMFBD) 

The speckle interferometry method presented above relies on a statistically average influ- 
ence of the wavefront aberration. In this Section we shortly present another approach that 
we have also used in this work. It is based on the simultaneous estimation of the object 
and the aberrations in a maximum likelihood sense using diff"erent simultaneous channels 
and several speckle frames. For more information see e.g. ( |Lofdahl|2002[|van Noort et aT 
20051 ILofdahl et al.|2007D . 



The method used is called Multi Object Multi Frame Blind Deconvolution (MOMFBD), 
which historically is a modification of the "Joint Phase Diverse Speckle" image restora- 
tion. The original method is based on the possibility of separating the aberrations from 
the object if we observe simultaneously in two channels introducing a known aberration, 
like defocussing the image, in one of them. Mathematically, both phase diversity and 
multi-object methods are particularizations from the "Multi Frame Blind Deconvolution". 
Using a model of the optics, including its unknown pupil image, it is possible to jointly 
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calculate the unaberrated object and the aberration, in a maximum likelihood sense. 



Coming back to Eq. 3.8 for a single isoplanatic speckle subfield, the Optical Transfer 
Function (OTF) is the Fourier transform of the Point Spread Function (PSF), which is the 
square modulus of the Fourier transform of the pupil function (P), that can be generalized 
with an expression like 

P = A ■ exp(i(f)) , (3.15) 

where A stands for the geometrical extent of the pupil (A= 1 inside pupil, A= outside). 
This unknown phase can be then parametrized using a polynomial expansion: 



= ^ O^mlAm, (3.16) 



meM 



where if/m, m e M, is a subset of a certain basis functions. The MOMFBD uses a com- 
bination of Zernike polynomials ( Noll||197 6 ) for tilt aberrations and Karhunen-Loeve for 
blurring effects, as they are optimal for atmospheric blurring effects ( Roddier[T990) . The 



{am} coefficients have therefore the information of the instantaneous wavefront aberra- 
tion, whether it comes from seeing conditions, telescope aberrations or AO influence. It 
is interesting to note that the expansion of the phase aberration is therefore finite (m e M) 



in our calculation, that leads to a systematic underestimation of the wings of the PSF (van 
Noort et al.|2005l ) 



For the calculation of the solution, the MOMFBD code uses a metric quantity that de- 
pends only on the {am} parameters and is expressed as the least square difference between 
the j speckle data frames, Dj, and the present estimated synthesized data frame, obtained 
by convolving the present estimation of PS F and object. 

J I y-' O* OTF P 

UM = V [ V ID/ - ^YJC "' ] (3.17) 

where the y term accounts for the noise and corresponds to an optimum low pass filter 
( Lofdahl|2002 ) and the u index for the spatial index in the Fourier domain. 



This mathematical expression, from Paxman et al. (1996), to solve the blind decon 



volution problem depends on the noise model used. In our case the MOMFBD assumes 
additive Gaussian statistics, which gives the simplest form of L and the fastest code, and 
turns to be appropriate for low contrast objects. 

The solution of the problem of image reconstruction is to find the set of {am} that min- 
imizes the metric L({Qr„,}), providing an estimation of the OTF, and from there the new 
estimation of the objects. Details on the process and optimization used can be found in 
Lofdahl ( 2002p . The final converging solution provides thus the real object and instanta- 



neous aberration simultaneously. 

With only one channel the {a^} are independent, but if we can specify linear equality 
constraints (LEG) to these parameters we can reduce the number of unknown coefficients 
for multiple channels. 

The Phase Diversity method is one example of LEG. By defocussing one of the cam- 
eras on a simultaneous channel we introduce a known phase contribution in the expansion 
of Eq. |3.16 This creates a set of related pairs of {a,,,}- Typically, 10 or even less realiza- 



tions of such pairs of images are enough for a good restoration. 
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Different channels observing simultaneously in different, yet close, wavelengths can 
be used also to constrain the {am}, as the instantaneous aberration can be considered the 
same for all channels. In our case we have several speckle images per position and two 
simultaneous channels. The broadband channel and the narrow-band channel scanning 
the spectral line at 21 positions with 20 frames per position. We have therefore a set of 21 
pairs of 2 simultaneous objects, with 20 frames for each object and channel. 

One interesting outcome of this multi object approach is that, if the observed data 
frames are previously aligned using a grid pattern, the resulting images are then perfectly 
aligned between simultaneous channels, which greatly reduces possible artifacts on de- 
rived quantities as Dopplergrams or magnetograms. 

In this work we have used this MOMFBD approach to process the data where our 
usual speckle interferometry method was not applicable. This mainly applies for on-limb 
observations, as the limb darkening gradient on the field of view influences the statistics. 
Also, with the actual presence of the off-limb sky, the data are not suitable for the narrow- 
band speckle reconstruction, as we don't have a broadband counterpart for the emission 
features present off the limb. 

The limb data set was reduced using this code (see Sec. 5.2), as well as some other 
data frames for comparison purposes with the speckle interferometry (Sec. 4.3 1. 

The MOMFBD code was implemented by van Noort et al. (2005 ) and was made freely 
available at www . momfbd . org. Given the high processing power needed it is written and 
greatly optimized in C++. It is developed to run in a multithreaded clustering environment, 
where the work is split in workunits and sent back from the slave machines to the master 
once the processing is done. A typical run with one of our Ha scans in broad and narrow- 
band channel, reconstructing the first 50 Karhunen-Loeve modes, takes ~ 7 hours to 
process with 20 CPU cores of 3.2 GHz. 



3.4 Infrared spectrometry 

For this work we have also used spectroscopic data in the infrared region, to study the 
spicular emission in the Hei 10830 A multiplet. For this purpose we used the echelle 
spectrograph of the VTT and the Tenerife Infrared Polarimeter (TIP). 

In this Section we summarize the instrument characteristics, the optical setup and the 
observations performed for the study of the emission profiles of spicules, which will be 
presented in Chapter [5] 



3.4.1 Instrument 

TIP was developed at the Instituto de Astroffsica de Canarias ( [Martmez Pillet et al.|1999| ) 
and recently upgraded with a new, larger infrared CCD detector ( jCoUados et al. 2007[ ). 
It is able to record simultaneously all four Stokes components with very high spectral 
resolution in the infrared region from Ifim to 2.3/zm, with a fast cadence and very high 
spatial resolution along the slit. 

The optical setup of the instrument is shown in Fig. 3.10[ After the main tank and 
the AO system, a narrow (~ 100 jum wide) slit is mounted in the plane of the prime focus 
of the telescope. The light reflected from the slit jaws enters a camera system to provide 
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Figure 3.10: Optical schema of the Tenerife Infrared Polarimeter (TIP) with slit jaw cam- 
era, predisperser and spectrograph of the VTT. After the AO correction, the light from 
the prime focus of the telescope enters the instrument through the slit. The light reflected 
from the slit jaws is recorded with video cameras to create context frames. After the slit, 
the polarimeter with the ferroelectric liquid crystals modulates the polarization of the light 
beam. The predisperser selects, with mask (pi), the spectral region to observe, and the 
spectrograph disperses the light into its spectral components. The nitrogen-cooled CCD 
detector records the modulated polarization of the spectra, dl and d2 are the diffraction 
gratings. 



images, to point the telescope and to have the region of interest imaged onto the slit. The 
small fraction of light entering the slit goes through the polarimeter, where the Stokes 
components are modulated. Then, the predisperser and spectrograph decompose the light 
into its spectral components. At the end of the optical path the detector is mounted, a 
CCD cooled below 100 K to reduce the thermal excitation of electrons in the CCD pixels. 

The polarimeter 

TIP is able to obtain simultaneously the full set of the four Stokes parameters that deter- 
mines the polarization of the light, from each point in the slit. However, this work concen- 
trates only on the intensity measurements. The polarization measurement is performed by 
means of two ferroelectric liquid crystals (FLC). These are electro-optic materials with 
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fixed optical retardation, whose axis can be switched between two orientations by apply- 
ing voltages of approximately + lOV. This amplitude of the rotation of the retardation axis 
is somewhat dependent on the temperature, and is ~ 45° at 20 - 25C. With two PLCs, 
with two possible states each, we can create four different combinations of modulation of 
the incident light. The four modulated intensities are four different linear combinations 
of {I,Q,U,V} with different weights on each parameter. With four consecutive measure- 
ments we can therefore retrieve the four components of the Stokes vector. Thus, TIP is 
able to obtain simultaneously the four components of the polarization for each full cycle 
of the polarimeter. Although TIP makes a full cycle of the PLCs in less than one second, 
we have to accumulate several spectrograms in order to increase the signal to noise ratio, 
especially when measuring weak signals like the polarization of spicules outside the solar 
limb. 

In the sequence following the light path, the physical setup of the polarimeter consists 
of a UV-blocking filter to protect the PLCs from intense high energy radiation at short 
wavelength. Then, the first PLC with a retardation of All and the second PLC with A/ A 
follow. The retardances of A/l and A/ A are nominal values. The actual retardances differ 
from these values and depend on wavelength. Pinally a Savart plate splits the light into 
two orthogonal linearly polarized beams. 

As part of the instruments we need a calibration optic subsystem (see explanation in 



Sec. 3.4.3) to account for the influence of the mirrors following the telescope. Por this 
reason, in front of the AO system, there is a polarization calibration unit (PCU) that can 
be moved into the light path. It is composed of a retarder with nominal retardance of A/ A 
in the optical spectral range, and a fixed linear polarizer. The retarder rotates a full cycle 
with measurements taken every 5 degrees, creating a set of 73 modulations of the light 
beam that are used to model the influence of the optics behind the telescope, but including 
AO, till the detector. The influence of the coelostat mirrors and the telescope proper on 
the polarization state are taken into account with a polarization model of these parts by 



Becketal. (2005). 



3.4.2 Observations 



Table 3.2 summarizes the details of the observing campaign for the course of this work. 
It focuses on studying the emission profiles observed in spicules in the Hei 10830 A 
multiplet. 

The strong darkening close to the solar limb and the presence of the limb make it 
difficult to use KAOS for off-limb observations, since the correlation algorithm of KAOS 
was not developed for this kind of observations. 

We scanned the full height of the spicule extension, starting inside the disc. We made 
a single spatial scan with long integration time per position. As the lock point of the AO 
was placed on a nearby facula inside the disc was chosen. Apart from the facula used for 
AO tracking, it was a quiet Sun region. In the present work we study only the intensity 
component of the Stokes vector (see definition in e.g. |Chandrasekhar| [T9601 |Wikipedia 
Stokes parameters]) ■ 
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Date 


Dec,4'^2005 


Location 


NE limb 


Spectral sampling # 


10.9 mA/px 


Time span 


1 scan in 66 min. 


Slit 


40"x or.'S 


Integration time 


5x2.5 s 


Step size 


a'35 


Max. height ofF-limb 


7" 


Seeing condition (ro) 


~ 7cm (max 12 cm) 


KAOS support 


yes 



Table 3.2: Characteristics of the data taken with TIP used in this work, ro is the Fried 
parameter. 



3.4.3 Data reduction 

As for the G-FPI case, the data reduction process aims to remove the instrumental effects 
as well as the atmospheric influence. For TIP data this involves three steps. The first 
is common to all CCD observations and consists in removing instrumental eff"ects, the 
second is the polarimetric calibration of the signal, and the third is the spectrosposcopic 
calibration. 



Reduction of CCD effects 

This processing is basically the same for all CCD observations: removal of dark counts 
and correction for diff"erential sensitivity of the pixel matrix with the gain table (using the 
flat fields). The only difl^erence to G-FPI data reduction is when creating the flat fields. 
The mean flat field frame is not flat. Although being a spatial average, it still contains 
spectral information. To retain only the gain table information we divide the flat field by 
the mean spectrogram, so that only the diff"erential response of the pixels is left (see Fig. 
3.1 ip . The mean spectrogram is obtained by averaging the flat field spectrograms over the 
spatial coordinate. 



Polarimetric calibration 

The signals recorded with the CCD are not directly the Stokes parameters (see description 
in e.g. Chandrasekhar| I960] ) . With two FLCs we have four different combinations in one 



full cycle. For each configuration in the cycle, we measure intensities as a particular linear 
combination of {I,Q,U,V} with different weights, so we can solve the ensuing system 
of equations. Also, in each CCD frame, we measure light of two orthogonal linearly 
polarized beams (see Sec. 3.4. 1[ ). 



An important problem in polarimetric observations is that each reflecting surface of 
the telescope changes the polarization state of the incoming light. So the optical path, with 
all the reflecting surfaces from the coelostat to the CCD, introduces a complex modulation 
of the incoming polarization. At the VTT there is a polarization calibration unit (PCU) 
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I 




(a) Frame of raw data 



(b) Flat field 




(c) Dark frame 



(d) Reduced frame 



Figure 3.11: Examples of the standard data reduction process for spectral data. The Flat 
field frame (b) is calculated dividing average flat field data by the mean spectra of the 
average. 

mounted in front of the AO system. This device feeds the subsequent optical components 
with light of well defined polarization states. So, once we have a set of Stokes parameters 
from different configurations of the PCU, we can obtain the modulation induced by the 
optical path, the Mueller matrix M, from the PCU to the polarimeter: 



/ 

Q 
U 
I V ) 



polarimeter 



I 

Q 
u 

y ) 



(3.18) 



input 



The inverse matrix of M will therefore relate the polarization state of the light that reaches 
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the polarimeter with the light arriving at the PCU position. However, the light path from 
the coelostat to the PCU (in front of the AO) cannot be calibrated with this system, so the 
reduction routines use a theoretical model of this part of the telescope. 

This process is already implemented with available reduction pipelines. Further in- 
vestigation of crosstalk or other additional polarimetric reduction are needed to reduce 
the instrumental effect in our data. However, this is not necessary for our case, since this 
work concentrates only on the intensity component. 



Spectroscopic reduction 

The last type of reduction procedure is related to the nature of spectroscopic data and 
consists of the calibration in wavelength, the continuum correction and a low pass filtering 
to remove noise. 

To calibrate our spectrograms in wavelength we make use of the two telluric lines 
present in our spectral range of the TIP data. Solar lines are subject to Doppler shifts from 
local flows and solar rotation. Yet, telluric absorption lines are formed in the atmosphere 
of the Earth. Therefore, they are always narrow due to only small Doppler broadening and 
are located at fixed wavelength. This provides a fixed reference coordinate that we use 
with the FTS atlas ( Neckel|1999[ ). Comparing both spectra we can accurately measure the 



spectral sampling which is for all data sets 10.9mA/pixel . See wavelength scale abscissa 
ofFig.|TT2j 



The transmission of the filters is not a constant in the transmitted wavelength range, 
so this creates an intensity variation curve in all our spectrograms. For normalization, we 
have to find the correct level of the continuum intensities of the spectrograms observed on 
the disc. For this, we use several spectral positions between spectral lines and calculate 
the ratio between the observed data and the values from the FTS atlas. We interpolate to 



create the continuum correction (see green dashed line on Fig. 3.12 1. 

An electronic signal was also found in some observed spectrograms with a frequency 
higher than those containing information on the solar spectrogram. We used for all data a 
low-pass filter which removes the power at all frequencies higher than a certain threshold, 
preserving the spectral line information. 

Once we have filtered and corrected the signal for all instrumental effects we have to 
remove finally the scattered light. We define the position of the solar limb as the height 
of the first scanning position (counting from inside the limb outwards), where the he- 
lium line appears in emission. For increasing distances to the solar limb a decreasing 
amount of sunlight is added to the signal by scattering in the Earth's atmosphere and 
by the telescope's optical surfaces. Since the true off-limb continuum must be close to 
zero, i.e. below our detection limit, the observed continuum signal measures the spurious 
light. Therefore, we removed the spurious continuum intensity level by using the infor- 
mation given by a nearby average disc spectrogram. This first subtraction estimates the 
continuum level on a region 6 A away from the Hei 10830 A emission lines. After this 
correction with a coarse estimate of the spurious light, a second correction was applied 
to remove the residual continuum level seen around the emission lines. This was needed 
since the transmission curve of the used prefilter is not flat but variable with wavelength. 



43 



3 Observations 




0828 



0830 



0832 
A (A) 



0834 



0836 



Figure 3.12: Example of intensity calibrated spectra on the disc near the limb. Raw 
spectrogram (blue line) has to be corrected for the continuum level to agree with the values 



in the FTS atlas ( iNeckel||1999[ black line). Using the continuum at several positions 
we can estimate the continuum correction (green dashed line). The corrected data (not 
filtered) are shown in orange. For the wavelength calibration we use the two telluric H2O 
lines (labeled in the figure). The region of the Hei 10830 A multiplet is also labeled, as 
well as some other lines in the range (Si, Cai , Nai). 
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chromosphere[] 



Since the discovery of the chromosphere 150 years ago, it has remained a lively and 
exciting field of research. Especially the chromosphere of active regions exhibits a wealth 
of dynamic interaction of the solar plasma with magnetic fields. The literature on the solar 
chromosphere, and on stellar chromospheres, is numerous. We thus restrict here citations 
to the monographs by |Bray and Loughhead] ( |1974[ ) and |Athay| ( |1976[ ) and to the more 
recent proceedings from the conferences Chromospheric and Coronal Magnetic Fields 
( |Innes et aL] [2005) and The Physics of Chromospheric Plasmas ( [Heinzel et al.||2007] ). 
With the latest technological advances we are able to scrutinize this atmospheric layer in 
great detail. The G-FPI in combination with post-processing techniques used in this work 
aims for the study of the temporal evolution of the chromospheric dynamics with very 
high spatial, spectral and temoral resolution. 

In this Chapter we present our investigations with the G-FPI inside the solar disc. The 
first Section focusses on data set "mosaic" and the presence of fast moving clouds. The 
subsequent Section presents the results of the investigation of fast events and waves from 
dataset "sigmoid". Finally we make a comparison between SI-l-AO and BD methods. 



4.1 Dark clouds 



As already noted in Sec. 1.2 the chromosphere is highly dynamic. Within and in the 



vicinity of active regions the interaction of the plasma with the strong magnetic fields 
gives rise to specially complex phenomena with fast flows. As an example we refer to a 
recent observation of fast downflows from the corona, observed in the XUV and in Ha 

Fast horizontal, apparent displacements of small bright blobs 
were observed in Hot by van Noort and Rouppe van 



by ITripathi etal] ( |2007] ). 
with velocities of up to 240 km s" 



derVoort(2006) 



Observations and data processing 



In this Section we use the data set "mosaic" (See Table 3.1 1 recorded on May, 31, 2004 
by K. G. Puschmann, M. Sanchez Cuberes and F. Kneer. It consists of a wide mosaiqued 
FoV around the active region AR0621. For each single FoV a series of five consecutive 
scans was performed, spanning a total of 4 min to study the temporal evolution. The FoV 



Contents from this Chapter have been partially published as Sanchez-Andrade Nuno et al. (2005 1, 



Sanchez- Andrade Nuno et al. ( 2007 1 
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of a single exposure was ~ 33"x23". To study a wide area the telescope was pointed con- 
secutively to 13 overlapping contiguous areas. The resulting mosaic covers a wide region 



with a total FoV of ~ 103"x94". In Figs. 4.1 and 4.2 we present the broadband image 



and narrow-band line core filtergram, respectively. In all mosaics, both in broadband and 
in all the narrow-band images there is a blank central area, that just corresponds to a small 
non-covered area. After dark subtraction and flat fielding, the data were processed using 



the SI approach (see Sec. 3.3.3). 



After the SI reconstruction, we have applied a destreching algorithm between the con- 
secutive broadband images to remove residual seeing eff"ects. The deformation matrix for 
the destreching was calculated for the broadband channel using a mean image as refer- 
ence. The same deformation matrix was then applied to the narrow-band spectrograms. 
To constrain the difi'erent frames of the mosaic of the broadband data, i.e. for connecting 
the individual subfields, a cross-correlation algorithm has been developed. The frames 
were smoothed by a boxcar of 5x5 pixels to take into account only large structures for 
the destreching and to reduce noise. The overlapping regions between the individual sub- 
fields have been used to scale the intensities and the several areas have been connected 
after proper apodisation. The arrangement of the individual subfields inside the broadband 
mosaic have been directly applied to the narrow-band data. 

Data analysis and interpretation 

We report the observations of numerous fast moving dark clouds in the FoV. Doppler- 
grams reveal that these clouds correspond to downward motion. Here we show a par- 
ticular fast dark cloud. Neither the continuum image nor the line center exhibit strong 
activity. However, if we study the filtergram taken in the red wing of the Ha line, a group 



of dark features becomes apparent (see panel 1 of Fig. |4.3[ ). 

Successive spectrograms every 45 s of the same region (panels 2 to 5 of Fig. |4.3| ) 
reveal a fast differential motion of this dark cloud. The position and direction is marked 
by the white arrow in Fig. 4.2[ A horizontal surface velocity of ~ 90 km/s is measured. 



Interestingly, the cloud has suddenly disappeared and was not longer seen in the last two 
observed frames. 

In Fig. |4.3| we display the corresponding spectral profiles for the central part of one of 
the cloud members (marked by white crosses in Fig. |4.3| ) at diff'erent times. 

We interpret the observed dark cloud, seen as a line depression in the red wing of 
the Ha line, as a signature of the Doppler shifts related to the fast movement of the 
dark cloud. From the spectral distance between the line core of Ha and the minimum 
position of the line depression we estimate a LOS downflow speed of ~ 51 km/s. This, in 
combination with the observed horizontal velocity leads to an approximate total speed of 
~ 103 km/s directed downwards. Further, the sudden disappearance of the cloud from the 
last 2 frames could be explained with a very strong related Doppler shift, thus the position 
of the line depression is displaced outside the scanned wavelength range. 
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Figure 4. 1 : Mosaic of speckle reconstructed broadband images of the active region NOAA 
AR0621, at yu = 0.68. The achieved high resolution by means of the adaptive optics and 
post factum reconstruction is ~ 0.2". The total area covered is ~ 103"x94". Limb is 
located to the left lower corner. 




Figure 4.2: Ha line center filtergram. It corresponds to one of the 1 8 reconstructed images 
along the spectral line. The resolution in these narrow-band images is < 0.5". One notes 
the various chromospheric features: ubiquitous short fibrils with different orientation, a 
wide bright plage region full of facular grains on the lower central part, and dark fibrils 
packed together outlining the magnetic field lines between sunspots around the central 



data gap. White arrow indicates position and direction of the dark cloud in Fig. 4.3 
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Figure 4.4: Broadband image of part of the active region AR10875 on April 26, 2006 at 
heliocentric angle § = 36°. The rectangles, denoted by A, B, B', C, and D, are the areas 
of interest (AOIs) to be analyzed and discussed below. 

4.2 Fast events and waves 

We continue investigating the active chromosphere on the disc of the Sun. We report 
on fast phenomena and waves observed in the Ha line with high spatial, temporal, and 
wavelength resolution. 



4.2.1 Observations and data reduction 



The observations correspond to dataset "sigmoid" in Table 3.1 They consist of a time 



sequence of 55 min duration of Ha scans with a mean cadence of 22 sec from the active 
region AR 10875. The observations were supported by the Kiepenheuer Adaptive Optics 
system (KAOS, von der Luhe et al.|2003[ ) under extremely good seeing conditions. 

Due to a technical problem, an increasing delay between successive scans was noticed 
during the observations. When the accumulative delay reached around seven seconds a 
new scanning procedure was restarted to avoid higher gaps between frames. This oper- 
ation needs around one minute. During the 55 minutes of this series, such an interrupt 
occurred twice, at 08:10:19 UT and 08:29:46 UT. This programming bug was corrected 
afterwards for future observations. 
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Figure 4.5: Narrow-band image corresponding to Fig. 4.4 in Hor at -I-0.5A offline center 



The same areas of interest are indicated as in Fig. 4.4 by the rectangles. 



The reduction process with SI+AO is explained in Sec. 3.3.3 We achieve a spatial 
resolution of ~(y.'25 for the broadband images at 630 nm and better than 0'.'5 for each of 
the 21 narrow-band filtergrams. Further, to follow the temporal evolution in time, both 
broadband and narrow-band images where cropped to the same common FoV, removing 
overall image shifts due to residual seeing effects. Afterwards, the speckle reconstructed 
broadband images were co-aligned to spatially and temporally smoothed images via a de- 
stretching code provided by Yi et al. ( 1992j ). The destretching matrix from the broadband 
image was also applied to the simultaneous narrow-band scan. To minimize the effects of 
the irregular sampling rate, the time sequences were interpolated to equidistant times with 
the cadence that leads to a minimum shift in time for each frame. This corresponds to a 
regular time step of 22 s. The data gaps at the times when observation was interrupted 
were filled by linear interpolation between closest observed images. 
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Figures 4.4 and 4.5 give the broadband scenery at t = 40.9 minutes during the series 
and the associated Ha image at +0.5 A offline center, respectively. The whole region was 



very active with a flare during the observation of the time sequence ( Sanchez- Andrade 



Nuno et aL]|2007b| ). The data set is certainly rich of information on the dynamics of the 



active chromosphere, especially since the spatial resolution is high throughout the se- 
quence. For the present study, we restrict further analyses and discussions to few regions. 
The areas of interest (AOIs) are indicated by rectangles and denoted by A, B, B', C, and 
D. In the presentations below the images from the AOIs were rotated to have their long 
sides parallel to the spatial co-ordinate in space-time images. AOI A contains a region 
where a long fibril developed twice during our observations. It has the appearance of a 
small surge ( Tandberg-Hanssen [1977 1. AOIs B and B' show a simultaneous fast event, 
possibly 'sympathetic' mini-flares with strong, small-scale brightenings in the Ha line 
core which last only few tens of seconds. AOIs C and D, with their long fibrils, are suit- 
able for the study of magnetoacoustic waves along magnetic field lines. Area C contains 
in its right part a region from which Ha fibrils stretch out to both sides and which, at the 
beginning of the time sequence, contained a small pore that disappeared in the course of 
the observations. Note also from Fig. 4.4 that the fibrils on the upper left side of area D 
originate in the penumbra of a small sunspot. 



4.2.2 Physical parameters 

The possibility to extract information from a good part of the Ha line profile in two 
dimensions and along the time series is highly valuable. We are interested in the physical 
parameters of the Ha structures. The line-of-sight velocities vlqs can be retrieved using 
the lambdameter method, while also many other parameters like the temperature and the 
mass density can be inferred by means of the cloud model. 

The lambdameter method ^Tsiropoula et al.|1993| ) is a common procedure to measure 
line of sight (LOS) velocities. It compares the Doppler shift of a spectral line with the 
position of the quiet Sun profile. We measure the profile bisector at several line widths. 
The method consists in measuring the displacement between the bisectors of the spectral 
profile and the reference quiet Sun profile. As pointed out by Alissandrakis et al. ( 1990| ) 



the resulting velocities give systematically lower LOS velocities than the cloud model (see 
below) by a factor of approximately 3. However, the qualitative behavior of both meth- 
ods are the same. The lambdameter method is therefore a fast method for a qualitative 
description of the velocity pattern of a region. 

The cloud model yields a non-LTE inversion technique. The formation of line profiles 
is the result of a complex interaction between the plasma and the radiation. Among oth- 
ers, the formation of a spectral line depends on the local temperature, velocity, chemical 
composition, magnetic fields, radiation fields, . . . Inversion techniques aim at retrieving 
this set of parameters from a given spectral profile. These techniques modify, in an it- 
eration scheme, the starting guesses of parameters based on certain assumptions until a 
converged solution, with modeled radiation and observations in close agreement, is ob- 
tained. We assume then that the calculated parameters producing the synthetic profile are 
the same as in the observed structure, as long as the assumptions are considered valid. 
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The cloud model allows the application of an in- 
version technique in cases when one can describe 
the radiation transfer through structures located 
high above the unperturbed solar photosphere. This 
method was first described by Beckers ( |1964[ ) and 



f to observer 



has been extensively used afterwards, e.g. by 



Tsiropoula and Schmieder 


( 1997]), Tsiropoula and 


Tziotziou|(2004), 


Tziotziou et al. (2004D, Al et al. 



(|2004V See also the recent review by Tziotziou 
( |2007 ). 



Figure 4.6 depicts the geometry of the cloud 
model. The considered "cloud" is located above 
the underlying photosphere at a height H and mov- 
ing at a speed t?. From the observer's position we 
can measure the projected proper motion relative to 
the background and the LOS velocity {Vios in Fig. 




||olar surface 



Figure 4.6: Geometry of the cloud 
model. 



4.6 1 as Doppler shifts. The observed intensity /(A/i) 
is the combination of the absorption of the back- 
ground intensity /()(A/l) with the emission from the 
cloud, dependent on the optical thickness of the cloud t(A/1): 



Jo 



(4.1) 



where S is the source function, which depends on the optical thickness along the cloud. 
In the model we make the following assumptions: 

1 . The structure is well above the underlying unperturbed chromosphere. 

2. Within the cloud, the source function 5, velocity and Doppler width are constant 
along the LOS. 

3. The background intensity profiles entering the cloud from below and in the sur- 
roundings are the same. 



to 



These assumptions simplify Eq. 4.1 

/(A/l) = Iq(AA) ■ e 
This, in terms of the contrast profile. 



+ 5(1- e-''-^'^) . 



CiAA) : = 



I{AA)-Io{AA) 
lo(AA) 



can be rewritten as 



C(AA) = i^—-l](l-e-'^^''). 



(4.2) 



(4.3) 



(4.4) 



,/o(Ai) 

Further, neglecting coUisional and radiative damping of the Ha absorption profile within 
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the cloud, the optical depth can be given by a Gaussian profile, i.e. 

( AA-AAj y 

t(AA) = Toe~y^) , (4.5) 

where tq is the line center optical thickness. Also, the central wavelength of the profile 
can be displaced due to a LOS velocity v of the cloud with a Doppler shift, AAj = Aqv/c, 
where Aq is the rest central wavelength and c is the speed of light. The width of the profile 
AAd depends on the temperature T and the microturbulent velocity trough the relation 



Ao 2kT 

AAn = -J (4.6) 

c y m 

where m is the atom rest mass. 

With these assumptions we end up with an inversion problem with 4 parameters: 5, 
AAd, Tq and vlos- ^^d is, in turn, the combination of 2 physical parameters (temperature 
and microturbulent velocity). 

More complex cloud models have recently been developed. These mainly focus on 
the nature of the source function S , allowing the parameters to vary along the LOS or 



multi-cloud models. However, as pointed out by [Alissandrakis et aL] ( |1990p , a simple 
Beckers cloud model like the one described above and used here provides useful, reason- 
able estimates for a large number of optically not too thick structures, tq < 1, for which 
the assumptions are adequate. 

In this work we also used the inversion in Ha structures where possible. The undis- 
turbed reference profile Io{A) is taken from a nearby area with low activity outside the 



FoV shown in Fig. 4.4 As the region under study was 'clouded out' in Ha, i.e. covered 
with structures to a large extent, the cloud model inversion failed often. In these latter 
cases, instead, the LOS velocity maps were determined with the lambdameter method or 
from difference images at Ha+0.5 A off line center with appropriate scaling. Calibra- 
tion curves to estimate from such Doppler-grams the true velocities were calculated by 



Georgakilas et al. ( 1990). From these we obtained that the velocities from the difference 
images were lower by a factor 2-4 than the true velocities, in agreement with those parts 
in the FoV where the cloud model inversion was successfully applied and with the results 



by |Tziotziou et al.| ( |2004l ) 



With the application of the cloud model and the inferred values of S , AAd, tq and 



Vlos we can derive other physical parameters. Following the approach by e.g. 'Tsiropoula 



and Schmieder ( 1997 ) we can calculate the population densities of the hydrogen levels 1, 
2, 3 (Ni,N2,Nj,), the total hydrogen density (Nh, including protons), the electron density 
(Ne), the total particle density (Nt) the electron temperature (Tg), gas pressure (pg), total 
column mass (M), mass density (p) and degree of ionization of hydrogen (xh)'- 



= '^^^^ (4.7) 

A^2= l.ieiQP'-^ cm"^ (4.8) 

A^3 = fM^ + ' (4.9) 

A^, = 3.2 10^VA^ cm"^ (4.10) 
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Nh= 5 10^ VA^ 

Nr= N, + il+ a)NH 

Pg = kNtT, 

M = (NnmH + O.OSSINh ■ 3.97m//)J 
p= Mid 
Xh = NJNh 



(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 



where d is the path length along the LOS through the structure, a is the abundance ratio 
of helium to hydrogen 0.0851), g2,g3 are the statistical weights of the hydrogen levels 
2 and 3 respectively, h is the Planck constant, v is the frequency of Ha, c the speed of 
light and k the Boltzmann constant. 



Table 4.1 summarizes average results from the cloud model and derived quantities for 



the long fibril in Fig. 4.7 



Parameter 


Av. value 


Parameter 


Av. value 


V [km/s] 


11.7 


Aiz) [A] 


0.34 


S/Ic 


0.154 


T 


1.05 


N2 [cm-3] 


4.5 • 104 


N, [cm-3] 


6.8 • 10'° 


[cm-3] 


1.1 • 10'' 


Ni [cm-3] 


3.8 • 10'° 


N3 [cm-3] 


4.2 • 10^ 




1.51 • 104 


p [dyn cm"^] 


0.38 


M [g cm-2] 


1.39-10-4 


p [g cm"^^] 


2.3 • 10"'^ 


Xh 


0.64 


Cs[km/s] 


14.4 







Table 4. 1 : Several derived parameters from the cloud model for the lower half section of 



the long fibril in Fig. 4.7 at ? = 25 min. We assume a LOS thickness equal to the width of 
the fibril (cylindrical shape) of 590 km and a micro-turbulent velocity of 10 km/s. First 
two rows result from the inversion technique while the others are parameters derived from 
them. 
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Figure 4.7: Space-time image of surge in AOI A at Ha + 0.5 A off line center, starting at 
14.7 min after the beginning of the sequence. The spatial axis runs along the minima of 
the surge intensities at this wavelength. 



4.2.3 Fast events in Ha 



Small recurrent surge 

Ejecta from low layers of active regions, called surges, have been observed in time se- 
quences of Ha filtergrams since many decades (e.g., |Tandberg-Hanssen|1977] ). 

In AOI A, a small surge occurred during the observed time series. It started near the 



pore at the upper right end of region A (cf. Figs. 4.4 and 4.5). It was straight and thin, 
with a projected length at its maximum extension of at least 15 Mm and with widths of 



approximately 2" at its mouth and 1" at its end. Figure 4.7 shows the temporal evolution 
of the surge in Ha -I-0.5A off line center. The space-time image starts 14.7 min after the 



beginning of the series and goes to the end of it. Along the spatial axis in Fig. 4.7, the 
minimum intensities along the surge are represented. 

The surge consisted of very thin fibrils, at the resolution limit < 0'.'5, being ejected 
in parallel. It started with several small elongated clouds lasting for 1-2 min. After- 
wards, it rose, reaching a projected length of around 14 Mm, and fell back after ~ 7min. 
Then it suddenly rose again after two min reaching lengths out of the FoV (more than 
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15 400 km) and lasted another five min before retreating again. And finally, the process 
recurred a third time, yet with lower amplitude in extension and velocity than for the 
first two times. The (projected) proper motion of the tip of the surge reaches a maxi- 
mum velocity of approximately 100 km s \ for both the ascent and the descent phases. 
Especially the second rise and fall showed large velocities. It is unlikely that the rapid 
rise and appearance of the surge in Ha are caused by cooling of coronal gas to chromo- 



spheric temperatures. The cooling times are much too long, of the order of hours (Hildner 



1974). Thus, the proper motions represent gas motions. The LOS velocities measured 



from Doppler-grams and corrected with the calibrations described above in Sect. 4.2.2 
amounted to +15 kms"^ during the ascent of the surge and reached -45 kms"' at the 
mouth during retreat. These latter velocities are lower than the proper motions. It thus 
appears that the chromospheric gas is ejected obliquely into the direction towards the 
limb. Average physical parameters in the surge obtained with the cloud model inversion 



are listed in Table 4.1 They are very similar to those of other chromospheric structures 



(see e.g. |Tsiropoula and Schmieder|1997[ ). 

Surges are known to show a strong tendency for recurrence, but on time scales of ~ 1 h. 



Sterling and HoUweg ( 1989|) have treated numerically rebound shocks in chromospheric 



fibrils and presented results in which a single impulse at the base of the involved magnetic 
flux tube drives a series of shocks on time scales of approximately 5 min. This appears to 
be a viable mechanism for the small surge observed here, apart from the initial conditions. 
The small 'firings' at the beginning of this surge suggest magnetic field dynamics that 
ultimately do cause a strong impulsive force, after some minor events. 



Synchronous flashes 

In the AOI pair (B, B') with a projected distance of ~14 Mm, brightenings occurred 
52.2 min after the start of the series in both sites at least as simultaneously as we can 
detect with the observational mode of scanning the Ha line. AOI B' is located in the 
umbra of a small spot with a complex penumbra and AOI B next to a pore. In between the 
two AOIs the sigmoidal filament ended while more structures of the extended and active 



filament system crossed the region between the two AOIs. Figure 4.8 shows the temporal 
evolution of the brightenings. 

The upper row of this figure is from AOI B', the lower from B. Two scans through 
the Ha profile are presented, of course without interpolation of the images to an identical 
time. The horizontal axes contain the run in both time and wavelength. 

The flash-like brightenings lasted only for less than 45 s, they were present neither 



in the scans before nor after the two scans shown in Fig. 4.8 The simultaneity of the 
two flashes, or mini-flares, suggests a relation between them. Possibly, one sees here a 
kind of sympathetic flares. These were discussed earlier in the context of synchronous 



flares excited by activated filaments ( Tandberg-Hanssenl|1977 ). Another interpretation is 



that one sees a mini-version of two-ribbon flares with a common excitation in the corona 
above them and simultaneous injection of electrons into the chromosphere. 

In AOI B, the flash exhibited sub-structure and apparently moved during the first pre- 
sented scan with speeds up to 200 kms"^. This strong brightening between 15 and 22 s 
has disappeared in the following scan. 

Figure |4^ depicts the recorded Ha profiles at the positions of the flash in AOI B, 
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as indicated by the arrows on the left side of Fig. |4.8[ The profiles are compared with 



those from the quiet Sun and from the average background. The profile from the isolated 
bright blob at 3.7 Mm (see inset in Fig. |4.9| ) shows a blue shifted emission above the 
background profile. This emission is still present in the following scan. At 2.8 Mm the 
line core is filled resulting in a contrast profile with strong emission (cf. Eq. 4.4). The 
profile at 1.1 Mm exhibits a strong emission beyond the continuum intensity in the red 
wing while the whole profile is enhanced above the background profile, the position of 
the emission peak would indicate a down flow with LOS velocity of 35 kms ^ It was 
shown by Al et al. ( 2004[ ) that such emission (contrast) profiles can be understood if one 
assumes an injection, likely from the corona, of much energy and electrons to obtain a 
response of the Hor line to temperature. These last two emissions at 2.8 Mm and 1.1 Mm 
have disappeared at the time of the following scan. 

Obviously, such fast events as in AOIs B and B' lie beyond the observing capabilities 
of our consecutive scanning method. We could however retrieve high spatial resolution 
filtergrams at several wavelengths to follow the temporal evolution at time scales of few 
seconds. With the present data set at our hand, we cannot decide whether the apparent 
proper motion of the flashing structure in AOI B is indeed as high as 200 km s"^ or whether 
the temporal resolution is too fast for the consecutive scanning. For example, the Hor 
profile from 1.1 Mm could have been in emission over the whole profile, but only for few 
seconds. It is, however, possible to design adequate observing sequences with duration of 
few seconds per scan, on the expense of taking filtergrams at fewer wavelength positions. 
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Figure 4.9: Ha profiles from the flash event. Each profile corresponds to an average 
over three pixels around the three selected points where the emission is highest in the 
blue wing, at central wavelength, and in the red wing respectively, corresponding to the 



white arrows in Fig. 4.8 (left half) at x=[3.7, 2.8, 1.1] Mm, respectively. For comparison 
the quiet Sun profile is also shown. Background profile corresponds to the mean from 
previous and following scans, where no brightening was found. The emission profile at 
1.1 Mm reaches an intensity of 1.1 of the quiet Sun continuum intensity. 



4.2.4 Magnetoacoustic waves 

In this Section, for the investigation waves in the chromosphere, we first refer to previ- 
ous observations of magnetoacoustic waves and then outline the magnetohydrodynamic 
(MHD) approximation (Ferraro and Plumpton 1966[ Kippenhahn and Mollenhoff| 1975[ 
Priest|1984) . From this, dispersion relations for sound waves, atmospheric waves, Alfven 
waves, and magnetosonic, or magnetoacustic, waves are derived by linearization. We 
continue discussing the observations of waves in long Ha fibrils. Finally we try an inter- 
pretation in the framework of waves in thin magnetic flux tubes. 



Previous observations of magnetoacoustic waves 

Apart from oscillations in sunspot umbrae (Beckers and Tallant 1969 Wittmann 1969| ) 



and running penumbral waves (e.g., von Uexkiill et al. 1983 and references therein). 



waves in the chromosphere were observed by many authors. E.g., Giovanelli ( 1975| ) 
describes waves along Ha mottles and fibrils with speeds of 70 kms~' and interprets 
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them as Alfven waves in magnetic flux tubes with approximately 10 Gauss field strength. 
Kukhianidze et al. ( 2006| ), from time sequences in Ha ofi" the limb, found kink waves 
in spicules with periods of 35-70 s. |Hansteen et al.| ([2006|), and |Rouppe van der Voort 



et al. ( 2007[ ) observed spicules and fibrils in the quiet Sun and in active regions with high 
spatial and temporal resolution observations. They succeeded via numerical simulations 
in explaining the dynamics of these chromospheric small-scale structures by magnetoa- 
coustic shocks, excited mainly by the solar 5-min oscillations (see also the simulations by 
Sterling and Hollweg|1989| and the review by |Carlsson and Hansteen|2005] ). 



Waves in the corona have as well been observed: E.g., Robbrecht et al. (2001 ) report 
on slow magnetoacoustic waves in coronal loops observed in high-cadence images from 
SoHO/EIT and TRACE. The speeds amount to 100 kms^^ |De Moortel et al.| (|2002l). 



also from high-cadence 171 A TRACE images, find that 3- and 5-min oscillations are 



common in coronal loops. They are also interpreted as magnetoacoustic waves. Tothova 
et al. ( 2007| ), from SoHO/SUMER data, study Doppler shift oscillations identified as slow 
mode standing waves in hot coronal loops. Fast-mode, transverse, incompressible Alfven 
waves, with speeds of 2 Mms \ in the solar corona were reported by Tomczyk et al 



(2007). 



Magnetohydrodynamic (MHD) approximation 

We use here the Gauss system of units. The MHD approximation is obtained from 
Maxwell's equations and the equation of mass conservation, the equation of motion, and 
an equation of state, under the following conditions: 

1 . The gas velocities v are small compared to the speed of light c, v «; c. 

2. Any changes are slow, such that phase velocities Vph = L/t «; c, where L is a typical 
length scale and t a typical time scale. 

3. The electrical conductivity cr is always very high such that the electric field is very 
small compared to the magnetic flux density, l^l \B\. 

4. One usually adopts in addition, to a good approximation, I) = ^ and § = ft, i.e. 
magnetic flux density and magnetic field have the same strength, in Gauss units. 

With j electrical current density. Maxwell's equations are then reduced to 

An ^ 

VxB=—j, (4.17) 

c 

VxE = --^, (4.18) 
c at 

V-i = 0. (4.19) 

Ohm's law, conservation of mass, and the equation of motion read as 

/=(r|£ + ^vxfij , (4.20) 
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^ = V-(pv^, (4.21) 

ot 

dv 1 
p^+p{v-V)^ = pi-Vp+-]xB. {AIT) 
ot c 

Here, p is the mass density, p the gas pressure, and g the gravitational acceleration vec- 
tor. The viscous, centrifugal, and Coriolis forces were omitted in the equation of motion, 
Eq. |4.22[ The equation of state relating the gas pressure p with mass density p and tem- 
perature T, is 

p = p{p,T). (4.23) 
Furthermore, we assume for simplicity adiabatic motion 

= 0, (4.24) 



dt \py 

with the ratio of specific heats 7 = 5/3 for monoatomic gases. 

The current density /can be eliminated by means of Ohm's law, Eq. 4.20 which yields 
the equation of motion 

P^+p(v-V)^/ = p^-V;? + -(Vx^)x^, (4.25) 
and the induction equation 



dt 



= -Vx|^Vx5j + Vx(vx5) , (4.26) 



where the last term on the rhs of Eq. 4.25 contains Maxwell's stress tensor. 

In an atmosphere with constant temperature and with the gravitational acceleration 
opposite to the vertical direction {g = -ge^, unity vector into z direction, \e^\ = 1) the 
hydrostatic equilibrium is 

-7- = -Pog , (4.27) 
dz 

with the solution 

Po = const! ■ ; p,, = const! ■ e'^'^ , (4.28) 
where the scale height is A = po/ipo ■ g). 



Magnetoacoustic gravity waves 

We consider now, to arrive at a dispersion relation for magnetoacoustic gravity waves, 
small perturbations from the equilibrium 

B = Bo + Bi; E = Eo + Er, /=/o + /i; (4.29) 

P = Po + Pi; p = Po + pi ; V = vo + vi . (4.30) 
Assuming further that Bq is homogeneous, vq = 0, and the conductivity is infinite. 
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cr ^ oo, one obtains 



_ ^ 47r ^ ^ ^ 

V X 5o = — jo = and (tEq = , i.e. Eq = . 

c 



(4.31) 



Inserting then the perturbed quantities from Eqs. |4.29] and [4.301 into the MHD equa- 
tions and neglecting quadratic terms and terms of higher order we obtain the linearised 
MHD equations 

+(;;j.V)po+Po(V-Vi) = 0, 



dt 



po- 



dp I 
dt 



+ {yi-V)p^-c] 



dp 
dt 



- + (vi-V)po 



= 0, 



dB\ „ /_, \ 
Po— = V X (vi X Bo) 



dt 



V-5, =0, 



(4.32) 
(4.33) 
(4.34) 

(4.35) 
(4.36) 



with the sound speed = [(7Po)/Po]''^^- From these one arrives, after some algebra, at 
the wave equation for the velocity 



dt^ 



c^V (V • i?i) - (r - 1) g gVvu 
+^[Vx{Vx(i?iX5o)}]5i/(4;r) . 



(4.37) 



Wave modes 

With Eq. |4.37| we make the ansatz 

v\{r, t) = vi exp ( A: • r - cufjj , 



(4.38) 



with wavevector k. 

For Bo = and ^ = one gets pure sound waves w ith phase velocity Vph = oj/k = c^. 
With ^0 = and ^ > one obtains atmospheric waves (Bray and Loughhe ad|l974 ). 

When the gas pressure is negligible, p = 0, and with g = 0, but |^ol > 0, the dispersion 
relation results 



<x»^vi/v^ = k^cos^avi - (k ■ vi)k cos a Bo + 



{k ■ Vi) - k cos a (Bo ■ Vi) 



k. (4.39) 



Here, Bo is a unity vector parallel to Bq, a is the angle between the wavevector k and Bo, 
and Va is the Alfven velocity with 



4;rp() 



, Pm,0 

' Po 



(4.40) 



where the magnetic pressure is Pm = B^/{Sn). Scalar multiplication of Eq. 4.39 with Bo 
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shows that Bo ■ vi = 0. This means that the (perturbed) velocity is perpendicular to Bq 
(since the Lorentz force on the pertur bed gas is perpendicular to Bq). 



Scalar multiplication of Eq. [4.39 



pertu] 
. 14J9 



with k yields 
-fcX)(^-Vi) = 0. 



(4.41) 



From this equation one can derive two magnetic wave modes: 



(1) Ass uming V • vi =0 gives the so-called incompressible mode and from the ansatz 
one gets k ■ Vi =0. Thus, this mode is a transversal mode with the 



4.38 



Eq. 

velocity perpendicular to the direction of propagation. From Eq. |4.39| we have 
£l) = icosavA. The waves are also called shear Alfven waves. For a = 0° one 
derives that Bi and vi are parallel and the propagation is along Bq. 



(2) Another solution of Eq. 4.41 h to = kvA, independent of a. These waves are com- 
pressional Alfven waves, and for a = 90° the velocity vi is parallel to k, i.e. we 
have longitudinal waves. 

Finally, admitting that the gas pressure is not negligible, p > 0, the phase velocity 
comes out as 



CO 



Vph = 



Uc-s^ + vl) ± (ct + v\- 2c,,vl cos 2a) 



1/2 



1/2 



(4.42) 



The '+' sign above gives the so-called 'fast magnetoacoustic waves' and the '-' sign the 
'slow magnetoacoustic waves'. Their phase speeds depend on a (see the hodographs in 
Ferraro and Plumpton|1966||Kippenhahn and Mollenhott11975||Priest|1984l ). 
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Figure 4.10: Average temporal power spectra of velocity from AOI C in Figs. 4.4 and 4.5 
before filtering (dashed) and after pass-band filtering (solid). 



Observational results 

The LOS velocities of the structures contain variations on long time scales of 10 min 
and longer as well as fluctuations with shorter time scales. To distill the latter, among 
them possibly magnetoacoustic waves, we applied a high-pass temporal filter and re- 
moved some high-frequency noise at the same time. The quantities then fluctuate about 
zero. Figure [4~T0l depicts the average power spectra of the LOS velocities in AOI C in 



Fig. 4.5 before and after filtering. We note that the 5-min oscillations are filtered out. 



while some oscillations at the acoustic cutoff" (corresponding to periods of approximately 



200 s) are partially retained. Yet the unfiltered and filtered power spectra in Fig. 4.10 do 
not show any predominant period. 



Figures 4.1 1 - 4.14 show examples of space-time slices from AOIs C and D. Tho fluc- 
tuations of several quantities are shown: 

1 . LOS velocities determined from differences of He intensities at +0.5 A off^ line cen- 
ter, henceforth referred to as Doppler-gram shces (bright indicates velocity towards 
observer); 

2. Ha line center intensities, henceforth LC slices; 



3. in Figs. 4.1 1 -4.13 diff'erences of intensities at -1-0.5 A offline center /o.5(?,+i)-/o.5(?,) 



with cadence At = tj+i - = 22 s, henceforth referred to as A/0.5 slices; 



4. in Fig. 4.14 diff^erences of intensities at line center Iiciti+i) - hciU), henceforth 
referred to as A/^c slices. 

Time runs from bottom to top with ? = at the start of the series. The interruptions/interpolations 
dXt ~ 18.0-19.6 min and 37.5-38.5 min are obvious. 
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Figure 4.11: Example of space-time slices, of I'.'l width, from AOI C in Figs. 4.4 and 
4.5[ From left to right: LOS velocity. Ha line center intensity, and intensity differences 
at Ha +0.5A off line center: /o.5(?,+i) - lo.si^i) with cadence of At = - = 22 s. The 
intensity differences in the right column are shifted up by 11 s. They are referred to as 
A7o.5 slices in the text. 
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Figure 4.14: Selected part of space-time slices from AOI D with slice width of 2'.'2. Left 
column: Ha line center intensity fluctuations; right column: intensity differences at Ha 
line center: ItciU+i) - hciU) with cadence of A? = ti+i - ti = 22 s, referred to as A/^c slices 
in the text. 



We focus attention to the oblique stripes in Figs .4.11-4.14 These are the signatures of 
magnetoacoustic waves. From their slopes we can measure phase velocities projected on 
the plane perpendicular to the LOS. In Fig. |4.1 Ij from AOI C, the waves appear to originate 
near the right edge of the AOI. This is one side at which the fibrils are rooted. Presumably, 
the waves are excited by the buffeting of motions at the photospheric foot points of the 
magnetic fields. As seen especially well in the Doppler-gram slices of Fig. 4.1 1[ but also 
in the LC slices, steep stripes originate from both sides of this region. The projected phase 
speeds are of the order of 8 km s ^ 

The stripes are often bent in the course of the temporal evolution, e.g. the wave 
parallel to the dashed line '1' in the A/0.5 slices of Fig. 4.1 1[ This wave starts off with a 
phase velocity of 14 kms"' and speeds up to approximately 40 kms \ one of the highest 
velocities measured. 

A prominent period is not detected. Sometimes, the waves appear to be repetitive, 
with two or three, at most, wave trains in sequence with periods between 90 s and 180 s. 
An example of consecutive wave trains is indicated by the three dashed lines '2' in the LC 
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slices of Fig. 4.1 1 Yet most time, the waves are solitary, with one single wave package 
traveling across the FoV. Many of the waves appear to spread out along the direction of 
propagation and to fade after having traveled a distance of 5-10 Mm. 

The amplitudes of the LOS velocities in the Doppler-gram slices are measured to 
approximately 1 kms~\ be it in the waves with low phase speeds or in those with high 
phase speeds. With the calibration discussed above in the context of the cloud model (see 
Sect. 4.2.2[ ) these amplitudes have to be multiplied with a factor of approximately 3. The 
resulting amplitudes are thus of the order of 3 kms"\ which is not a small perturbation 
compared with the sound speed (c.f. below the discussion on the magnetoacoustic waves). 

Figure |4.12| from AOI D shows similar space-time slices as those from AOI C in 
Fig. 4.11 Yet here, the waves are excited at both sides and travel into the AOI, sometimes 
crossing from left and right and possibly colliding as in the example parallel to the dashed 
lines '1' in the Doppler-gram. The long lasting (more than 7 min), solitary wave train 
(parallel to dashed line '2' in Fig. 4.12 1 has a phase velocity of approximately 13 km s~\ 
a typical speed of the 'slow' waves in this AOI. A correction for foreshortening, i.e. that 
we see only the projection of the phase speed on the plane perpendicular to the LOS is 
to be excluded since the fibrils in AOI D, as well as those in AOI C, are oriented almost 
perpendicularly to the direction to the limb (see Fig. |4.5[ ), thus perpendicularly to the 
LOS. The wave parallel to the dashed line '3' in the A/0.5 slices of Fig. 4.12 gives a phase 
speed of 30 km s again not to be corrected for foreshortening. This is a typical phase 
speed of the fast waves. 

Figure [4~T3] gives another example of space-time slices from AOI D, with wider slices 
of 2'.'2 width and shorter time span of 13 min duration than in Figs. |4. 11] and |4.12| The 



long lasting wave train in the Doppler-gram slices (parallel to dashed line '1') gives again 
the typical phase speed of 13.3 kms ' with (calibrated) LOS velocity amplitudes of ap- 
proximately 2 kms ^ These LOS velocities are transversal, in the sense that they are 
perpendicular to the propagation and to the Ha fibrils. The 'fast' wave in the A/0.5 slices 
(parallel to dashed line '2') exhibits also the typical phase speed of 32 kms~^ with cali- 
brated LOS velocities of approximately 1.5 kms^^ 

Figure [4T4l give s a 7.25 min long section of the temporal development of fluctuations 
in AOI C with slice widths of 2'.'2, but this time the LC slices and the A/^c slices only. 
Note that dark and bright features in Ha LC indicate increased and decreased absorption. 



respectively, not enhanced and reduced temperature (see Al et al. 2004, van Noort and 



Rouppe van der Voort 2006). The two solitary waves between the pairs of horizontal 



dashed lines (a, a') and (b, b') have phase speeds of approximately 25 km s"^ Inspection 
of the LC slices shows that the waves consist of elongated, thin blobs with length of 1"- 
2" and width of approximately 0'.'5. Apparently, the waves do not travel in the spatial 
direction along straight lines, but along sinuous lines with deviations from straight lines 
of approximately 0'.'5 in amplitude. This suggests that on these small scales the magnetic 
field is not straight and homogeneous but entangled. 

The presentation of the difference slices A/^c in Fig. 4.14| is prepared to study temporal 
displacements of absorption features. These are only seen if the displacements have a 
strong component perpendicular to the LOS. The bright and dark small-scale features, 
lying parallel and next to each other, in the upper part of the A/^c slices, between the 
dashed line pair (b, b') are suggestive of such displacements perpendicular to the direction 
of propagation. 
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We summarize in short the observational findings on magnetoacoustic waves: 

1. Generally, we find two kinds of waves: slow waves with phase velocities of 12- 
14 km s"' and fast waves with phase velocities of 25-33 km s"' (maximum velocity 
found 42 km s The waves appear to develop from low phase speed to high phase 
speed waves and vanish after having traveled a distance of 5-10 Mm. 

2. h-respectively of the wave mode, the LOS gas velocities are of the order of 2- 
4kms~^. 

3. The waves are mainly solitary. They consist of short (l"-2") and thin Of.' 5) blobs 
of compressed gas. 

4. The waves appear to follow wiggly, entangled magnetic field lines with possible 
lateral displacements. 



Interpretation - waves in thin magnetic flux tubes 



For the interpretation of the observations from AOI C and D, we adopt the picture of 
waves in thin magnetic flux tubes, whose radius is small compared to the pressure scale 
height. The propagation of waves in magnetic flux tubes were treated by, among others. 



Defouwl ( fT976l ), |Wentzel| ([1979]), |Spruit| (T9ST ), and recently by [Musielak et al] dlOOT] ). 
Spruit assumes a thin, cylindrical magnetic flux tube parallel to the z axis, with radius 
R, magnetic field along the tube of strength B, gas pressure p, mass density p, and tem- 
perature T. The gravity is neglected. The tube is embedded in an external medium with 
properties B^, p^, p^, and T^. Inside and outside the tube the magnetic and atmospheric 
parameters are constant. In Spruit's ( 1982[ ) work, the MHD equations are linearized and 
a mode analysis is performed, with proper conditions at the interface between flux tube 
and surrounding medium. 

Incompressible Alfven waves (V • vi = 0, with small velocity perturbation vi) are also 
possible in flux tubes. They are torsional Alfven waves. The compressive solutions lead 
to 

V • vi = A S,„(nr) exp[z {ojt + mcp + kz)] , (4.43) 

with amplitude A, !Bm(nr) Bessel functions of order m, r the distance from the axis of 
the tube, and (p the azimuthal angle. Inside the tube, the waves propagate along the z 
direction. For n the relation holds 



= (co^ - vlk^) {0? - c]k^)IV{a? - c^) (v^ + c^)] 
Here, the tube speed Cf is introduced with 



2 2 



(4.44) 



(4.45) 



which shows that the tube speed is smaller than both the Alfven and the sound velocity. 

Spruit] ( 1982) showed that in the limit kR the mode with m = is a longitudinal 
mode with Vph = Ct, which is approximately the sound speed c.v for va » c^. This mode 
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is often referred to as the 'sausage mode', with velocity inside the tube parallel to the 
magnetic field. 

In the same limit and for m > one obtains the so-called 'kink waves', with phase 
speeds related to the magnetic fields and densities through 



P+Pe 



1 B^ + B^ 
An p+ Pe 



(4.46) 



These waves are transversal waves, and Spruit's (19821) analysis takes into account the 
dragging by the ambient medium. The phase speeds are obviously v^^^ = 



VI for Pe = 



p. 



5, 



5; Vj;, 



v^/2 forp^ = p, Be = 0, and v^;, = 2 • forp^ = 0, B 



B. 



We now compare the observations of waves with the expectation from this linear wave 
theory. We adopt that the waves propagate along the magnetic field and that the influence 
of gravity on the wave properties is negligible. The period at the acoustic cutofl^ of 200 s 
is longer than the periods, actually seen only rarely, in our data. Likewise, the period for 
the cutoff of kink waves ( Spruit|1981[|Choudhuri et aL]|1993 1 is approximately 400 s, for 
small plasma /3, which is the ratio of gas pressure to magnetic pressure, /3 = (Snp)/B^. 

With the parameters in Table |4.1| for the surge discussed above in Sect. |4.2.3| i.e. with 
gas pressure p = 0.38 dyncm"^ and mass density p = 2.3 x 10 gcm"^, the sound 
velocity is = 16.6 kms"^ From the determination of parameters in a wide range of 
chromospheric Ha structures by [Tsiropoula and Schmieder| ( |1997| ), |Tsiropouia| ( |2000| ), 



and [Tsiropoula and Tziotziou] ( |2004[ ) we obtain values of the sound speed in the range 
of 13.5-16.7 kms~^ The widely found temperatures of T = 10^ K and the mean molar 



mass of 0.8 from the ionization equilibrium of hydrogen found from Table 4.1 and from 
the above works give a sound velocity of 14.4 kms The phase velocities of the slow 
waves observed here are compatible with these values, if one accounts for possible small 



projection effects and for a small reduction for the velocity of tube waves (cf. Eq. |4.45[ ). 
De Pontieu et al.| ( [2004| ) adopted magnetic field strengths of the order of 100 Gauss 



in the chromosphere of active regions. With this value and the commonly found mass 
densities of 0.8-2.3x10"^^ g cm~^, the Alfven velocity is = 1 000 . . . 600 km s"\ much 
higher than the velocities of the fast waves in the present observations. We believe, that 
100 Gauss is an upper limit of the field strengths in the chromosphere of AOIs C and 



D. From high spatial resolution (approximately O.'SS) data from a plage region by Bello 



Gonzalez and Kneer| (2008) we find an average field strength in the photosphere of 60- 
90 Gauss. This may possibly be reduced by a factor of 2 in chromospheric fibrils as in 
AOIs C and D by spreading out of the field lines over areas which possess little field in 
the photosphere. Otherwise the fibrils would not be so elongated. Yet still this yields to 
Alfven velocities of Va ~ 200 km s~\ as a minimum value. 

Giovanelli (1975) has measured velocities of 70 kms^ in chromospheric Ha struc- 



tures. With a magnetic field strength of 10 Gauss and with reasonable particle densities 
he arrived at the Alfven velocity in agreement with these measured phase velocities. In 
the present work, one would need field strengths as low as 5 Gauss for an Alfven velocity 
of 32 kms~^ as observed. We note that even with 5 Gauss the motions are still dominated 
by the magnetic field, i.e. «; 1 holds. 

We estimate the maximum phase speed measurable from our data to 250-300 km s~^ 
Such velocities would still be detectable. The phase speeds found here are in the range 
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25-35 kms~^. (The highest measured speed amounts to 42 kms"^). These are obviously 
incompatible with Alfven waves in a homogenous magnetic field with 30-100 Gauss. We 
mention several possibilities to reconcile our measurements with the picture of fast mode 
magnetoacoustic waves along the magnetic field, i.e. of Alfven waves. 

1. The magnetic field strength in the fibrils of AOIs C and D is indeed as low as 
5 Gauss, which is not very probable considering the very high activity in the whole 
area observed. AOIs C and D are not especially located at the outskirts of this 
activity. 

2. Propagation of a fast mode wave in a flux tube surrounded by a medium with low 
or zero field strength but with high gas density would reduce the phase speed (cf. 
Eg. [446]). 



Apparently, the waves start as slow mode waves with phase velocities of the order 
of 10-14 kms"^ and then are transformed into fast mode waves propagating with 
Alfven velocity. Yet the transformation does not occur immediately. Examples 
are seen in Fig. 4.1 1[ While the solitary waves evolve into fast mode waves their 



wave packages get dispersed and they decay by spreading out along the direction 
of propagation. 

We do not measure phase velocities but group velocities of solitary wave pack- 
ages. We have calculated for the slow mode the cusped surface of the wave front 



according to Ferraro and Plumpton (1966 cf. their Fig. 13), which is rotation 



ally symmetric about the direction of the magnetic field. The adopted Alfven and 
sound velocities were 200 kms ' and 16 kms~\ respectively. The maximum ve- 
locity of this surface is only marginally larger than the sound speed by 3.2%, and 
the maximum deviation from the direction of the magnetic field is (f.Ol. Thus, the 
propagation of such slow mode pulses is practically along the magnetic field with 
the sound velocity. 

The picture is actually more complicated: The waves with low phase speed seen 
here are not pure longitudinal waves. The gas velocities of the waves have a strong 
transversal component of the order of 3 km s^. Furthermore, the propagation of the 
fast waves deviates from straight lines, their motion appears more wiggly, possibly 
because the magnetic fields are entangled. Under the aspect of these observations 
the linear theory of small perturbations of straight flux tubes appears to be not suf- 
ficient. 
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4.2.5 Summary on some properties of the active chromosphere 

Thanks to the good resolution, we could follow the evolution of small-scale chromo- 
spheric structures of an active region. From the rich dynamical processes in the ob- 
served, very active, flaring region some areas were selected for detailed investigation in 
the present work: 

1. A small surge: It showed repetitive occurrence with a rate of some 10 minutes. 
The surge developed from initial small active fibrils to a straight, thin stucture of 
approximately 15 Mm length, then retreated back to its mouth to reappear again 
two times. The gas velocities reach approximately 100 kms"^. The rebound shock 
model by Sterling and Hollwegl ( 1989j) seems to be a viable explanation. 



2. Two small-scale, synchronous, possibly sympathetic flashes, or mini-flares: In a 
pair of small areas, two brightenings occurred simultaneously and disappeared dur- 
ing two Ha scans with total duration of 45 s. Presumably, the evolutionary time 
scale is much shorter, few to 10 s. Yet we could follow the evolution with a tempo- 
ral resolution of 2 s by analysing Ha filtergrams at difl^erent wavelengths. One of 
the two flashes showed an apparent proper motion with a speed up to 200 km s"^ 
while it was developping a high emission, above the continuum intensity, in the red 
part of the Ha profile. However, the cadence of the scanning was too slow to decide 
whether the temporal evolution consisted in a rapid horizontal proper motion with 
a final fast down flow or in a rapid change of emission at fixed local postitions. 

3. Magnetoacoustic waves in long fibrils: In two areas with long fibrils, the structures 
exhibited many magnetoacoustic waves running parallel to the fibrils, thus presum- 
ably also parallel to the magnetic field. The waves are mostly solitary. Few times, 
two or three repetitive wave trains could be seen with periods of 100-180 s. The 
waves start at the footings of the fibrils with a speed of 12-14 km s~\ which is not 
much lower than the sound speed estimated for such structures and similar to the 
tube speed. Most of the waves get accelerated to reach phase speeds of approxi- 
mately 30 km s ^ Then they spread out along the fibrils and fade. The final phase 
speed is much lower than the Alfven speed of > 200 km s^^ estimated from reason- 
able magnetic field strengths in the active region chromosphere of 30-100 Gauss 
and reasonable mass densities in the fibrils of 2x10"^'' g cm"^. Furthermore, we ob- 
serve that the slow waves have strong transversal (LOS) velocity components with 
~3 kms"\ i.e. are not purely longitudinal, and that the fast waves consist of short 
(l"-2"), thin (~0'.'5) blobs and apparently move along sinuous lines. We conlude 
from these findings that a linear theory of wave propagation in straight magnetic 
flux tubes is not sufficient. 
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4.3 Comparison between speckle interferometry and blind 
deconvolution 



In Sec. 3.1 we introduced the image degradation problem due to atmospheric distor- 



tions for all ground based solar observations and astronomical observations on general. 



In Sec. 3.2.1 we explained the adaptive optics approach used at the VTT to reduce the 



image distortions in real time. Finally, in Sec. 3.3.3 we described the basis of two differ- 
ent post factum image reconstruction techniques, the Speckle Imaging (SI) and one type 
of Blind Deconvolution (BD) with simultaneous Multiple Objects and Multiple Frames 
(MOMFBD). 

In the case of data from the G-FPI, the SI approach reconstructs separately the broad- 
band images and uses both the original recorded frames (after dark subtraction and flat 
fielding) and the reconstructed image from this channel to obtain the reconstructed narrow- 
band images at the various positions along the spectral line. 

The MOMFBD code applied to our G-FPI data uses at the same time, for each spectral 
position, the various (15) pairs of simultaneously recorded broadband and narrow-band 
frames. Thus, at 21 wavelength positions along the spectral line, two different objects 
were observed and their images were reconstructed with 15 frames per object. 

In this thesis work we used the SI method for fields of view on the solar disc. The 
last version of the code takes into account the field dependence of the PSF around the AO 
lockpoint, so we will refer to this version as SI-l-AO. However, data frames on and off the 
limb cannot be reconstructed with the current code. Near the limb the contrast is lower 
than near disc center, which makes any reconstruction more diflicult. Moreover, KAOS 
can lock on the low-contrast structures near the limb only under very good seeing condi- 
tions. Also, the limb darkening at large heliocentric angles makes it difficult to determine 
the STF on the rings around the lockpoint for the SI-i-AO. Beyond the problems inside 
the limb, off-the-limb emission features seen in the narrow-band images lack broadband 
counterparts, and therefore there exist no simultaneous data from which we could apply 
the second part of the SI to reconstruct the off-hmb parts of the images. 



Dataset limb (results in Sec. 5.2) was recorded under very good seeing conditions, 
with KAOS locked on a nearby facula correcting 27 (Zernike polynomial) modes most of 
the time. For the post factum image reconstruction we used the MOMFDB, for which the 
limb darkening presents no problems. Spicules above the limb do not posses a simulta- 
neously observed broadband object, so it is expected that their spatial resolution is lower, 
since there are no multiple objects, just the multiple frames for each spectral position in 
the narrow-band channel. 

In order to compare results from both approaches we have reconstructed with both 
methods the same field of view on the disc, a frame of the sigmoid data set. In this 
Section we present the comparison. As we show, for the BD case we made two different 
reconstructions with different limits of the expansion of the aberration in Karhunen-Loeve 
(K-L) modes. In the first case the expansion of the wavefront aberration is done until de 
n"^ K-L polynomial, while on the second run we expanded to the first 100 modes. The 
running time of the code is highly sensitive to the number of modes, being slower with 
more modes. 
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Figure 4.15: Results for different image reconstruction techniques. Axes are in arcsec- 
onds. White rectangles enclose the region where the RMS contrast is calculated. Upper 
left: Best speckle raw image, contrast is 5.9%. Upper right: SI+AO result, contrast 
is 11.1%. Lower left: MOMFBD running 17 modes, contrast is 6.7%. Lower right: 
MOMFBD running 100 modes, contrast is 10.0%. 



broadband comparison 



speckle frome 
burst average 
SI+AO 
MOMFBD, 7 
MOMFBD,, 
SI - M0MF3D, 




54.65 01.46 00.74 00.49 

anqulo^ scale (arcsec) 



00.37 



00.29 



Figure 4.16: Comparison of the power spectra (in arbitrary units) of the same image for 
different image reconstruction techniques. 
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Figure 4.17: Results for different image reconstruction techniques for the line center 
narrow-band filtergram. Scales on the axes are in arcseconds. Upper left: Best speckle 
raw image. Upper right: SI+AO result. Lower left: MOMFBD running 17 modes. Lower 
right: Image difference. In this case the differences reach 45% of the fluctuations in the 
reconstructed frames. 
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Figure 4.18: Comparison of the power spectra of the same image for different image 
reconstruction techniques. 
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Broadband 



Figure 4. 15 compares the full FoV image, while Fig. 4. 16 shows the corresponding power 
spectra. The speckle frame corresponds to the one with highest rms contrast (5.9% inside 
the white rectangle). The SI+AO reconstruction shows a much higher resolution, with 
more power at all frequencies than the speckle frame for angular scales larger than ~ 0'.'32. 
The reconstructed image has also less noise than the speckle frame (small-scale end of the 
power spectra) . The rms contrast of granulation for this image is 11.1%. In the case of the 
MOMFBD with 17 modes the power of the reconstructed image is significantly lower than 
for SI+AO, albeit having a lower noise level (comparable even with the burst average). 
We have to run up to 100 modes to arrive at a similar contrast as for SI+AO. The noise 
threshold for the last run, coincides with that of the SI+AO approach. The rms contrast 
of granulation for these images are 6.7% with 17 modes and 10.0 % with 100 modes. 



Figure 4.16 shows also the power spectrum of the difference between SI+AO and 
MOMFBD (green line), which is many orders magnitude lower than one of the power 
spectra themselves. Only at scales smaller than ~ 0'.'4, the difference becomes of the same 
order as the power spectra. Taking the differences of the reconstructed images shows that 
99.8% of the pixels in the FoV have intensity fluctuations lower than 15% of the intensities 
in the images themselves. 



Narrow-band 



The narrow-band images have lower intensities than the broadband images, especially 
at the core of the Ha line. Much less images are used for reconstruction, so a lower 



resolution is expected. Figure 4.17 compares the images at the Ha line center, while 



Fig. 4.18 shows the corresponding power spectra. The SI + AO reconstruction shows a 
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Figure 4.19: Close-up subfield from the narrow-band spectrogram at the Ha core. Axes 
are in arcseconds. 
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higher resolution than the speckle frame, with more power at scales larger than ~ Of! 5. 
At smaller scales, the speckle frame is dominated by noise. The MOMFBD with only 
17 modes gives already a similar resolution than the SI+AO and better treatment of the 
noise. Fig 4.19| shows a close-up region where the better noise treatment of the MOMFBD 
is clearly visible. 

The difference between the methods is bigger than in the broadband case, as expected 
since the intensity and resolution are lower. Nonetheless the agreement is very high, 
92.8% of the pixels in the difference image have amplitudes smaller than 0.15 of the 
average intensity in the reconstructed images (similar results are found for other spectral 
position, reaching 99% in the wings, at wavelengths +1 A off the line center). 



Conclusions 

In this Section we have shown the good convergence of both post-processing approaches. 
Using different techniques we arrive at similar results and spectral profiles. The ampli- 
tudes of the difference images are lower than 0.15 of the average image intensity in more 
than 99% of the broadband and above ~ 90% for the narrow-band images. In the case 
of the broadband reconstruction it was necessary to use 100 modes for the MOMFBD 
method to reach similar results as for SI-i-AO, while, in the narrow-band case, already 
with 17 modes the MOMFBD gives better images than SI-l-AO. 

The main disadvantage of BD methods is the computational load. The reconstruction 
of the single data set from broad and narrow -band and only 17 modes takes ~ 7 hours 
to process with 20 CPU cores of 3.2GHz. For the 100 modes run, given the limited 
resources, we only used the broadband frames (Multi Frame BD). If the data set "sigmoid" 
were reconstructed with BD methods, even with only 17 modes, it would have taken 
around 130 days on our computing resources. 

The main advantage of the BD is its ability to reconstruct an image even with only 
few frames. This is of special importance when observing fast evolving targets. The 
SI needs much more frames. The property of reconstructing simultaneously recorded 
images from different "objects" (e.g. broadband and the Ha narrow-band) leads to a 
perfect sub-alignment of the results, which avoids spurious signals in derived quantities. 
Note however that, not simultaneously observed objects, like in the several consecutive 
scans with the G-FPI, are not aligned since they are not recorded under identical seeing 
conditions. 

The SI-l-AO method is considerably much faster, around 10 and 15 minutes for the 
broadband and narrow-band images, respectively, with the same computers used for the 
MOMFBD reconstruction and gives better results for the broadband reconstruction, even 
using 100 modes in the latter method. However, with MOMFBD, the resolution and 
treatment of the noise is better in the narrow-band case. The main current advantage of 
the BD methods for our work and data is the possibility of reducing narrow-band limb 
and off-the-limb data scans. 

Anisoplanatism is an issue common to both approaches. In both cases the large FoV 
is divided into smaller subfields, where the assumption of isoplanatism is valid. It is 
therefore important to address this point for both cases. The image difference does not 
show any sub field pattern. However, there can still be some small effects. For this reason 
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we have used the integrated contrast profile of the diff"erence, defined as 

\Isi+AoW — iMOMFBDnW 



A 



iSI+AO 



(4.47) 



where Is i+aoW and Imomfbd„(^) correspond to the reconstructed images using the SI+AO 
method, and to the images using MOMFBD with 17 modes, respectively. 

CI qualitatively measures the total difference between the profiles. If they were equal, 
then CI would be 0, while an increasing difference in the profiles increases the value of CI. 
Since the subfield locations are the same for all the spectral positions, this information is 
added along the scan, while the intensities of the structures at each position are essentially 
subtracted out. The subfield pattern does not disappear with the subtraction of images 
reconstructed with different methods since they do not necessarily coincide. 



Figure 4.20 shows the calculated CI. The weak subfield pattern is revealed, especially 
in regions where the difference is low (dark background). The mean edge length of the 
squares is approximately around 32 pixels. 

The amplitude of the grid pattern is very low, only revealed after the calculation of 
CI. Presumably this comes from the apodization. When joining common regions on 
overlapping subfields, the common parts are overlaid in the final image. This, while 
preserving the structures, reduces the noise, which leads to slightly smaller noise levels in 
these overlapping lanes. This grid is common for all wavelength positions. The difference 
between the methods is low enough to reveal this small decrease of the noise level (leading 



to darker areas in Fig. 4.20) when the total effect is calculated, by using the CI parameter. 
Therefore, regions with more contrast, where the difference between SI and BD is bigger, 
the presence of this pattern is masked, as shown in the figure. 
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5 Spicules at the limb [] 



Spicules, known for more than 130 years (see the hand drawings by Secchi|1877 ), repre- 
sent a prominent example of the dynamic chromosphere. We refer the reader to reviews 
by |Beckersl ( |1968[[1972| ) and to the paper by |Wilhelm| ( |2000[ ) on UV properties. Accord- 
ing to these works, spicules are seen at and outside the limb of the Sun as thin, elongated 
features. They develop speeds, measured from both proper motion and Doppler shifts, 
of 10-30 km s"' and reach heights of 5-9 Mm on average, during their lifetimes of 3-15 
minutes. Recent observations from HINODE (Kosugi et al. 2007 p and own results pre- 
sented below in Sec. |5.2| have changed the traditional picture. Some spicules live for only 
few seconds, and spicules may be much more inclined with respect to the vertical than 
adopted hitherto. 

As pointed out by Sterling] (2000 1, a key impediment to develop a satisfactory under- 
standing has been the lack of reliable observational data. Many theoretical models have 
been proposed to understand the nature of spicules, using a wide variety of motion triggers 
and driving mechanisms. In this Chapter we focus on the Hei 10830 A triplet emission 
line (see Sec. |2.3[), using recent technical improvements in observational facilities, and 



on the results from the limb observations in Ha. 



5.1 Spicule emission profiles observed in He 1 10830 A 



The energy levels that take part in the Hei 10830 A triplet are basically populated via an 
ionization-recombination process ( Avrett et al.||1994| ). The EUV coronal irradiation (CI) 
at wavelengths A < 504 A ionizes the neutral helium, and subsequent recombinations of 
singly ionized helium with free electrons lead to an overpopulation of all ortho-helium 
levels. Alternative theories suggest other mechanisms that might also contribute to the 
formation of the helium lines relying on the coUisional excitation of the electrons in re- 
gions with higher temperature (e.g., [Andretta and Jones|[T997] ). We are able to provide 
observational evidence of the link between the corona and the infrared emission of this 
line, in the frame of the current theoretical models of the solar atmosphere. 



Centeno (2006) modelled the ionisation and recombination processes using various 
amounts of CI, non-LTE radiative transfer, and different atmospheric models (see also 
Centeno et al.|2007 ). They have simulated limb observations for different heights, obtain- 
ing synthetic emission profiles in spherically symmetric models of the solar atmosphere. 
One conclusion of their study is that the ratio of intensities CR = /wue/^red) of the 'blue' to 
the 'red' components of the Hei 10830 A emission is a very good candidate for diagnos- 
ing the CI. The population of the metastable level depends on optical thickness, whose 



Contents from this Chapter have been partially published as Sanchez-Andrade Nuiio et al. (2007a I 
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Figure 5.1: Measured Hei 10830 A emission profiles for increasing distances to the solar 
limb, scanning a broad range of the height extension of the spicules. Each profile is the 
average of the 312 pixels along the slit (which was always kept parallel to the limb). 



variation with height governs the change in the ratio "Rasa function of the distance to the 
limb. 



Trujillo Bueno et al. (2005) measured the four Stokes parameters of quiet-Sun chro- 



mospheric spicules and could show evidence of the Hanle effect by the action of inclined 
magnetic fields with an average strength of the order of 10 G. They modelled the Hei 
10830 A profiles assuming the medium along the integrated line of sight as a slab of con- 



stant properties and with its optical thickness as a free parameter. Trujillo Bueno et al. 



(2005 1 showed that the observed intensity profiles and their ensuing 91 values can be re- 
produced by choosing an optical thickness significantly larger than unity. |Centeno| ( |2006| ) 
demonstrated that this optical thickness is related to the coronal irradiance (through the 
ratio "R), thus providing a physical meaning to the free parameter in the slab model (see 
also Centeno et al. 2007). 



5.1.1 Observational intensity profiles and intensity ratio 

We present novel observations showing the spectral emission of Hei 10830 A and its 
dependence on the height of the spicules above a quiet region. We compare the deduced 
observational K with that obtained from detailed non-LTE numerical calculations using 
available atmospheric profiles. 



These data correspond to the data set described in Table 3.2 After the standard re 
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5.1 Spicule emission profiles observed in He i 10830 A 




Figure 5.2: 3D representation of the measured He 1 10830 A emission profiles for increas- 
ing distances to the solar visible limb. Note that the x-axis is wavelength, the y-axis the 
height above solar limb and the z-axis the intensity normalised to the maximum emission 
in the line center of the red component. 



duction process (Sec. 3.4.3 ) we obtain 21 intensity profiles above the infrared limb, with 



a step size of 0'.'35. Figure 5.1 shows the emission profiles of the Hei 10830 A (after 



the reduction process) for different heights above the limb. Figure 5.2 illustrates this in 
three dimensions, as a function of wavelength and the distance to the solar limb, clearly 
showing a change in the intensity ratio of the blue and red components of the multiple! 
(^ = /blue //red) with height. 

For the calculation of K we need to determine the amplitudes of the blue and red 
components of the emission profile (as shown in Fig. 5.3[ ). 

To determine the core wavelength of the red component of the triplet we fitted a Gaus- 
sian profile to its core, in a 1.3 A range around the maximum. After symmetrising the 
observed profile around this maximum, using the values on the red side of the red compo- 
nent, we fitted another Gaussian function to the resulting symmetric profile. Subtraction 
of the fitted symmetric profile from the data leaves the emission profile of the blue com- 
ponent, which was also approximated by a Gaussian to determine its central wavelength. 
Our tests trying to fit directly both profiles using two Gaussians failed in a number of 
cases, probably due to the following reasons: (a) the red component is in fact the result of 
two blended lines, (b) the much weaker blue component was almost hidden in the broad- 
ened red component, and (c) the presence of noise. Our technique determines first the red 
component and then, after subtraction of the fitted profile, the blue one. 
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Figure 5.3: Determination of the blue and red components of the Hei 10830 A triplet 
from the observed emission profiles. In this example the slit was placed at Y.'4 off" the 
solar visible limb. See text for details. The solid line represents the average emission 
profile. The dotted line is the Gaussian fit to the symmetrised red component. Subtraction 
of this from the observed profile leaves the blue component, which is also fitted by a 
Gaussian profile (thin solid line). The sum of both Gaussians (dashed line) gives the fit to 
the observed profile. 



We have thus separated the helium emissions into their red and blue components as- 
suming only that both are present and that they are both symmetric. We can now measure 
their widths and intensities and also check that the line core positions coincide with the 
theoretical ones. After the fitting process the residuals between measured and observed 
profiles were small, the largest errors occurring in the determination of the core inten- 
sities of the red line. This happens because the red component consists of two blended 
lines (with a separation of 0.09 A), a fact that flattens the emission profile near the core as 
opposed to a more peaked Gaussian function. Nevertheless, the diff"erences between fitted 
profiles and data are only significant in the red core and are always lower than -1-0.08 of 
the maximum normalised intensity, with a mean difference of ~0.02. To avoid system- 
atic errors, we used the observational values for the center of the red component when 
calculating !R. 



5.1.2 Results 

The chromospheric temperature and density are too low to populate the ortho-helium lev- 
els via collisions ( |Avrett et al.||1994p . The EUV irradiation from the corona (CI) ionises 
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the para-helium, and the subsequent recombinations lead to an overpopulation of all the 



ortho-helium levels, in particular those involved in the 10830 A transitions. Centeno 



(2006) and Centeno et al. (2007) have modelled the off"-the-limb emission profiles and 
concluded that the ratio 'R = hiuelhed is a function of the height and a direct tracer of the 
amount of CI. Here we compare the results from the theoretical modelling with observa- 



tions. 



Trujillo Bueno et al. (2005 ) modelled their spectropolarimetric observations assuming 
a slab with constant physical properties with a given optical thickness. In the optically 
thin regime 7? ~ 0.1 2, which is the ratio of the relative oscillator strengths of the triplet. 
As the optical thickness (at the line-center of the red blended component) grows, this ratio 
also increases until it reaches a saturation value slightly larger than 1 for r ~ 10. (This 
type of calculation can be done and improved as explained in Trujillo Bueno & Asensio 



Ramos 2007). To reproduce the observed emission profile Trujillo Bueno et al. (2005) 
had to choose r ~ 3. Interestingly, the values of r yielded by this modelling strategy are 
consistent with the more realistic approach of Centeno (2006), where non-LTE radiative 
transfer calculations in semi-empirical models of the solar atmosphere are presented, us- 
ing spherical geometry and taking into account the ionising coronal irradiation. With our 
data we are able to test such theoretical calculations by comparing the measured values of 
H with those resulting from various chromospheric models. This way we may eventually 
trace the amount of CI inciding on the spicules. The analysis described above yielded the 
values of K for the observed profiles. The resulting dependence on the distance to the 
solar limb, for each pixel along the slit and each position of the slit above the limb, are 



presented in Fig. 5.4 The solid black line gives the average value of H. 

The dependence of 'R with height can be understood in a qualitative way as follows: In 
the outer layers of the chromosphere the density is so low that the transitions occur in the 
optically thin regime. With decreasing altitude the ratio K increases (proportionally with 
density) until a maximum optical thickness is reached. At even lower layers, although the 
density still continues to rise, the extinction of the coronal irradiance leads to a reduction 
in the number of ionizations, which results in a decrease of the optical thickness in the 
core wavelength of the red component, and thus in a decrease of K. 

For a quantitative comparison with theoretical modelling we have used the results 



from Centeno (2006) and Centeno et al. (2007) where they calculated the ratios % for 



different standard model atmospheres: FAL-C and FAL-P ( Fontenla et al.|199T] ) and FAL 



X ( Avrett| 1995). The FAL-C and FAL-X models may be considered as illustrative of the 



thermal conditions in the quiet Sun, while the FAL-P model of a plage region. The FAL-X 
model has a relatively cool atmosphere in order to explain the molecular CO absorption 
at 4.6 yu. 



The comparison is shown in Fig. 5.5 We notice that the modelled height variations 
of K agree only in a qualitative manner with what is found in our observations. However, 
the calculations from difi'erent models of the solar atmosphere are unable to reproduce the 
measured ratio. Higher values of the coronal irradiance lead to an increase of the optical 
thickness (at the line centers of the He i multiplet) and an upward shift in the run of K vs. 
height. Yet the shape of the height dependence is mainly given by the atmospheric density 
profile and the attenuation of the ionising radiation as it reaches the lower layers of the 



chromosphere. It is also clear from Fig. 5.5 that the models do not extend high enough 
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Figure 5.4: Measured ratio 7? = hiuelh-ed as a function of distance to the solar limb. Thin 
lines come from each pixel along the slit. The thick solid line represents the average and 
the dashed line the value of the optically thin regime. 



5.1.3 Conclusions 

The theoretical behaviour of the ratio K agrees qualitatively with observations. Yet, a 
quantitative comparison shows poor agreement. Also, the simulated ratios are highly 
model dependent. As already explained, the failure to reproduce the observed profiles is 
very hkely due to the density stratification not being adequate for spicule modelling and 
to the limited vertical extension of the atmospheric models. Modelling of the intensity 
ratio in the He i infrared triplet should account for the fact that the solar chromosphere 
is inhomogeneous on small scales and that the spicules are small-scale intrusions of chro- 
mospheric matter into the hot corona. 



5.2 High resolution imaging of spicules 



De Pontieu et al.| (|2007') recently published high resolution observations of spicules with 



the Solar Optical Telescope on board Hinode ( |Kosugi et al.||2007] ) in the Call H line 
at 3968 A. They find at least two types of spicules that dominate the structure of the 
magnetic solar chromosphere: Type I with 3-7 minute timescales that correspond to the 
hitherto known spicules, and the new Type II spicules, developing in ~ 10 s, and lasting 
10-150 s. These are also very thin, with widths down to the spatial resolution (120 km). 
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Figure 5.5: Observed (average) vs. theoretical variation of the ratio 'R= hiuelhed with 
height. 



Also, Trujillo Bueno et al. (2005 1 used spectropolarimetric observations and a theoretical 
modeling accounting for radiative transfer effects. They find that the magnetic field in 
spicules is aligned with the visible structures and measure field strengths of up to 40 G 
with an inclination of 35° with respect to he local vertical. 

This Section compares these observational and theoretical properties with our high 
spatial resolution observations with the G-FPI. We use the dataset "limb". It consists 
of several consecutive Ha scans with a field of view that includes the limb and a region 
outside the disc up to a height where no emission from spicules in the Ha core is observed. 
The seeing conditions were extremely good during the observations and the AO system 
could lock on a nearby facula. After usual dark subtraction and flat fielding we have used 
the BD method (see Sec. 3.3.3.C I to achieve highest spatial resolution. We were observing 
the limb near both poles. Li Figures |5]6] and |5]7] we present some examples of the reduced 
data. 

Our time series of four minutes duration already shows a wide range of dynamics. We 
observe both long lasting spicules and fast evolving phenomena. Measuring the inclina- 
tion of the projected spicules to the local vertical we find angles up to 30° for the north 
pole, as it has been reported ( ,Beckers|1968 ) . The projected height above the limb varies 
from the wings to the core, from 2770 to 3750 km at +0.5 and +1A respectively. Near 
the south pole we find, however, much stronger emission and higher inclinations. The 
maximum angle is close to 70° from the local vertical, while the maximum height reaches 
up to 8250 km. We also find one horizontal fibril/spicule, as well as the presence of kinks 
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Figure 5.6: Reconstructed narrow-band scan observed near the solar north pole. The size 
of each image is 56'.' 1 x 19'.'!. The wavelength of the filtergrams increases by 0.1 A from 
top left to bottom right row by row. The third image in the third row is closest to the 
center of the mean line profile. The images have been rotated to have the limb parallel to 
the horizontal axis. 



or bends in some spicules. The width of single resolved spicules varies from a maximum 
width of 1 000 km at the spicule footing to a minimum size of 250 km, almost down to the 
resolution limit of the images, both in faint spicules and in others with strong emission. 
We also can retrieve the spectral profile at each pixel. 



Figure 5.6 demonstrates an important contribution to the understanding of spicules. It 



solves the long-standing question about the counterparts of spicules on the disc (|Grossman i 



Doerth and Schmidt] 1992 1. The first and last four filtergrams of the scan across Ha in this 
figure show that spicules outside the limb continue as dark fibrils inside the disc. 



In Fig. 5.8 we show the mean variation from the disc to the limb of the intensity in 



Ha around the north pole. Further, Fig. 5.9 presents mean intensity variations from the 
disc to the limb for several wavelengths around the Ha line center. The emission at the 
line center is almost constant from the disc up to a height of around 5" above the limb, 
where the intensity starts to decrease. 
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5 Spicules at the limb 




2 4 

arcsec obove vis. limb 



I 



Figure 5.8: Image representation of the mean measured spicular profiles from image 5.6 



the X-axis is the height above limb, while the y-axis is the wavelength around Ha line 
center (black horizontal line). Horizontal cuts at A - Aq = [0, +0.5, +1]A are shown in 
Fig. 15:9 




Figure 5.9: Average over lY.'l of Ha intensity profiles inside and outside the limb, for 
several line positions, observed near the solar north pole. Dotted line: broadband intensity 
at 6300 A, the inflection point defines the position of the solar limb; green thick line: 
intensity at Ha line center, which is nearly constant till «5" above the limb and then 
decreases outwards; dashed lines: intensities at -0.5 A (blue) and +0.5 A (red) off line 
center, with height of spicular visibility decreasing at « 3'.'7; solid lines: intensities at 
+ 1A off line center. It is seen that the Ha line turns from absorption inside the limb into 
an emission line (line intensities higher than the continuum intensity) above the limb. 
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6 Conclusions and outlook 



We have studied the dynamics of the solar chromosphere, both on the disc and above 
the limb, using two spectral regions (Ha in visible light and the infrared Hei 10830 A 
multiplet). By means of real-time correction and different post-processing techniques 
we have reduced the image degradation induced by the Earth's atmosphere achieving 
resolutions in Hor up to O'.'S and better. This Chapter summarize the main conclusions of 
this work. 



Observations and analysis 

The basic results from the observations taken from the disc are: 

• Data taken in the combination with the "Gottingen"-Fabry Perot Interferometer (G- 
FPI), adaptive optics and speckle interferometry have high quality. We have ob- 
tained a time sequence in Ha of 55 min from the active region AR 10875 at he- 
liocentric angle « 36°. The time cadence is 22 seconds, and its field of view 
77" X 94". For each interpolated time step we can retrieve 23 filtergrams along the 
Ha spectral line with 45 mA FWHM and spatial resolution better than 0'.'5 . Simul- 
taneous broadband images at 6300 A were also obtained, with spatial resolution of 
0'.'25, close to the telescopic diffraction limit. 

• We have observed the dynamics of a small surge in detail: It showed repetitive oc- 
currence with a rate of some 10 min. The surge developed from initial small active 
fibrils to a straight, thin structure of approximately 15 Mm length, then retreated 
back to its mouth to reappear again two times. The gas velocities reach approx- 



imately 100 kms The rebound shock model by Sterhng and HoUweg (1989) 
seems to be a viable explanation. 

• The region was very active during the observations. We studied two small-scale, 
synchronous, possibly related flashes, or mini-flares. The simultaneity is within 
seconds, while their total evolution time was ~ 45s. The brightenings were sepa- 
rated by ~ 14 Mm. The used scanning parameters of the G-FPI were slow for this 
fast evolution, yet we could follow it with a temporal resolution of 2 s by analysing 
filtergrams taken consecutively at different wavelengths across the Ha line. One of 
the two flashes showed an apparent proper motion with a speed up to 200 kms"\ 
while developing a high emission in Ha, above the continuum intensity. 

• For the observations of waves we restricted our study to two areas exhibiting long 
fibrils. Yet the results likely represent the typical behavior of chromospheric mag- 
netoacoustic waves within this active region. By means of high-pass frequency 
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filtering, we observe waves running parallel to the fibrils, thus presumably also par- 
allel to the magnetic field. They were mostly solitary waves, although sometimes 
repetitive wave trains could be seen with periods of 100-180 s. Most pulses start 
with velocities on the order of 12-14 kms~' and get accelerated to reach phase 
speeds of approximately 30 kms"^ Furthermore, we observe that the slow waves 
have strong transversal (LOS) velocity components with ~3 kms~\ i.e. are not 
purely longitudinal, and that the fast waves consist of short (l"-2"), thin (~0'.'5) 
blobs and apparently move along sinuous lines. 

Further, we have analyzed observations of spicules inside the disc and above the limb 
with the G-FPI data. Given the properties for this kind of observations we could not use 
the speckle interferometry method to reduce the atmospheric distortions. Instead we have 
used the blind deconvolution approach, in particular the version developed at the Swedish 
Institute for Solar Physics for multiple simultaneous objects with multiple frames per 
object ( |van Noort et al.||2005| ). The observations and analysis yielded the following main 
results: 

• It is possible to successfully use multi-object multi-frame blind deconvolution meth- 
ods with the G-FPI to reduce atmospheric distortions. This is specially important 
for on-limb observations, where the current speckle interferometry method is not 
applicable. 

• We have observed spicules in Ha at both solar polar caps. Compared with the solar 
north pole, we find much stronger spicular emission at the south pole that could be 
related to the presence of a coronal hole. The maximum projected height reaches 
8250 km, while we see inclinations of the spicules up to 70° form the local vertical. 
We can resolve the detailed structure of the spicules as well as the presence of kinks 
or bends in some cases. The width of a single spicule ranges from 1 000 km down 
to the resolution limit of around 250 km. 

• Using the retrieved spectral profile we can observe that spicules outside the limb 
continue as dark fibrils inside the disc, as shown in Fig. |5.6[ This answers a long 



standing question, e.g. cited by Grossmann-Doerth and Schmidt (1992). 



Thus, we have used two different post-processing approaches to reduce the image degra- 
dation with the Hot spectral line data. Since both methods are based on different ap- 
proaches, we have reduced the same observational data with both techniques and com- 
pared the results. These are the main conclusions: 

• The agreement of the images from both approaches is high. The achieved resolu- 
tion comes close to the diffraction limit in both cases. Even though both methods 
split the image into isoplanatic subfields for individual reconstruction, there is no 
difference of subfield re-composition when comparing the results. 

• In general, the biggest advantage of speckle interferometry over blind deconvolu- 
tion is the small computational time required. A complete restoration of one full 
Hot scan like the ones used in this work is roughly ~ 400 times faster with speckle 
interferometry than with blind deconvolution methods. 
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• The main advantage of blind deconvolution methods is their versatility. It can be 
applied with only few frames, with one or few simultaneous objects, or both at the 
same time. This method is highly advisable when aiming for e.g., fast evolving 
targets, or limb observations. The perfect sub-alignment of simultaneous objects 
avoids spurious signals on deduced quantities, like magnetograms. 

• For the broadband channel we find that the speckle interferometry gives images with 
high contrast. Only when forcing the blind deconvolution method to reconstruct up 
to high (100) Karhunen-Loeve modes, the results are similar. We note that the speed 
of the reconstruction is proportional to the number of modes. 

• The narrow-band images are clearly better reconstructed with the blind decon- 
volution, even with only 17 Karhunen-Loeve modes. The noise treatment gives 
smoother images with details at smaller scales (of the order of « Qf.'S). 

Further, we have obtained and analysed spectroscopic measurements in the infrared. We 
have centered our studies here on the intensity profiles of the Hei 10830 A multiplet 
above the limb. 



Recent work, like e.g. by Trujillo Bueno et al. (2005), Centeno (2006), has demon 



strated the importance of the intensity ratio between the blue and red component 
of this triplet as tracer of the coronal irradiance. In this work we present novel ob- 
servations showing the variation of this parameter with distance to the limb with a 



resolution of 0'.'35 up to 7" above the solar visible limb (See Fig. 5.4). 



Interpretation of observations 

For the interpretation of the observed data we have used several models and previous 
theoretical results to compare with the presented data. The main results from this analysis 
are: 

• From the intensity profiles of the Ha spectral line inside the disc we can infer many 
physical parameters. We have applied the lambdameter method as a fast and easy 



way to retrieve qualitative velocity maps. Also we have used Beckers's 1964| cloud 
model. Our simple non-LTE inversion code provides the possibility to infer many 
physical parameters, e.g. hydrogen and electron density, mass density, tempera- 
ture, . . . The results are in agreement with the data given in the current literature. 

• From the linearization of the MHD problem, we discuss the interpretation of the 
observed waves as magnetoacoustic waves. We assume estimates with reasonable 
magnetic field strengths in the chromosphere of the active region of 30-100 Gauss 
and reasonable mass densities in the fibrils of 2x10"'^ gcm~^. The observed wave 
speeds are much lower than the expected Alfven velocities. We conclude from these 
findings that a linear theory of wave propagation in straight magnetic flux tubes is 
not sufficient. 

• From the infrared observations we have calculated the ratio of amplitudes in the 

has modelled 



two main components of the Hei 10830 A multiplet. Centeno (2006 



synthetic limb observations according to the current theories of formation of this 
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triplet and chromospheric models. The agreement is only qualitative. The failure 
to reproduce the observed profiles is very likely due to the density stratification 
not being adequate for spicule modelling used and to the limited vertical extension 
of the atmospheric models. Modelling of the intensity ratio should account for 
the fact that the solar chromosphere is inhomogeneous on small scales and that 
the spicules are small-scale intrusions of chromospheric matter into the hot corona. 
Future models of the solar chromosphere should be constrained by the observational 
evidences presented within this work. 



Outlook 

The solar chromosphere represents a lively and exciting field of research. The wealth of 
structures, its dynamics and the wide range of evolution timescales are the consequence 
of the peculiar properties of this atmospheric layer. Current instruments like the ones 
used here, are able to observe and study in great detail new phenomena, that test current 
models and, as a consequence, helps our understanding of the solar atmosphere. This 
thesis aimed at contributing to the understanding. Yet, work to extend this research is 
already in progress. Here we shortly describe some of this work and give an outlook to 
further steps to be undertaken next. 

• The blind convolution method provides a practical way to study the spicules in Hor 
near and above the limb. Data from a short time sequence, taken under very good 
seeing conditions, are currently under reconstruction with phase diversity methods. 
The analysis will shed light onto the dynamic phenomena in spicules. 

• We have learned that the sequential scanning, with the G-FPI, with cadence of 22 s 
is not fast enough in some cases. For future observations, we can design scanning 
modes of 2-3 second resolution taking images at fewer wavelength positions in a 
spectral line, like Ha. 

• New infrared data of spicules near the solar poles and the equator, below coronal 
holes or in coronal active regions will help us to understand the detailed formation 
process of the He i 10830 A lines. 

• Full Stokes spectropolarimetric data of the Hei 10830 A multiplet are available 
from an earlier observing campaign. Scans above the hmb were performed under 
very good seeing conditions. We can therefore extend our present analysis and 



study the polarization. We aim to investigate the Hanle effect as suggested by Tru 



jillo Bueno et al. (|2005 ) and make use of available inversion techniques like e.g. 



from |Lagg et al.| ( |2004j . 



The new Gregor telescope ( [Balthasar et al.|2007| ) will host the G-FPI from the com 



ing year on. The combination of this new facility with other instruments like Hinode 
( |Kosugi et al.|2007| ), will provide new exciting resources for further research. 



94 



Bibliography 



Al, N., Bendlin, C, Hirzberger, J., Kneer, R, Trujillo Bueno, J., 2004, Dynamics of an 
enhanced network region observed in Ha, A&A, 418, 1 131-1 139 

Alissandrakis, C. E., Tsiropoula, G., Mein, P., 1990, Physical parameters of solar H-alpha 
absorption features derived with the cloud model, A&A, 230, 200-212 

Andretta, V., Jones, H. P., 1997, On the Role of the Solar Corona and Transition Region 
in the Excitation of the Spectrum of Neutral Helium, ApJ, 489, 375-394 

Athay, R. G., 1976, The solar chromosphere and corona: Quiet Sun, Reidel, Dordrecht, 

xn 

Avrett, E. H., 1995, Two-component modeling of the solar IR CO lines, in Infrared Tools 
for Solar Astrophysics: What's Next?, 15th NSO Sac Peak Workshop, (Eds.) J. Kuhn, 
M. Penn,pp. 303-311 

Avrett, E. H., Fontenla, J. M., Loeser, R., 1994, Formation of the solar 10830 A line, in 
Infrared Solar Physics, lAU Symp. No. 154, (Eds.) D. Rabin, J. Jefferies, C. Lindsey, 
pp. 35^7 

Balthasar, H., von der Liihe, O., Kneer, R, Staude, J., Volkmer, R., Berkefeld, T, Caligari, 
P., CoUados, M., Halbgewachs, C, Heidecke, R, Hofmann, A., Klvaiia, M., Nicklas, 
H., Popow, E., Puschmann, K., Schmidt, W, Sobotka, M., Soltau, D., Strassmeier, 
K., Wittmann, A., 2007, GREGOR: the New German Solar Telescope, in The Physics 
of Chromospheric Plasmas, (Eds.) P. Heinzel, 1. Dorotovic, R. J. Rutten, vol. 368 of 
Astronomical Society of the Pacific Conference Series, pp. 605-610 

Beck, C, Schlichenmaier, R., CoUados, M., Bellot Rubio, L., Kentischer, T, 2005, A po- 
larization model for the German Vacuum Tower Telescope from in situ and laboratory 
measurements, A&A, 443, 1047-1053 

Beckers, J. M., 1964, A study of the fine structures in the solar chromosphere, Ph.D. 
thesis. University of Utrecht 

Beckers, J. M., 1968, Solar Spicules (Invited Review Paper), Sol. Phys., 3, 367^33 

Beckers, J. M., 1972, Solar Spicules, ARA&A, 10, 73-100 

Beckers, J. M., Tallant, P. E., 1969, Chromospheric inhomogeneities in sunspot umbrae, 
Sol. Phys., 7, 351-365 



95 



Bibliography 



Bello Gonzalez, N., Kneer, R, 2008, Narrow-band full Stokes polarimetry of small struc- 
tures on the Sun with speckle methods, A&A, in press 

Bendlin, C, 1993, Hochauflosende zweidimensionale Spektroskopie der solaren Granu- 
lation mit einem Fabry-Perot-Interferometer, Ph.D. thesis, Gottingen 

Bendlin, C, Volkmer, R., 1995, The two-dimensional spectrometer in the German Vac- 
uum Tower Telescope/Tenerife. From observations to results., A&AS, 1 12, 371-382 

Bendlin, C, Volkmer, R., Kneer, P., 1992, A new instrument for high resolution, two- 
dimensional solar spectroscopy, A&A, 257, 817-823 

Berkefeld, T, 2007, Solar adaptive optics, in Modern solar facilities - advanced solar 
science, (Eds.) F. Kneer, K. G. Puschmann, A. D. Wittmann, pp. 107-114, Univer- 
sitatsverlag Gottingen 

Born, M., Wolf, E., 1999, Principles of Optics, Cambridge University Press, 7th edition 

Bray, R. J., Loughhead, R. E., 1974, The Solar Chromosphere, Chapman and Hall, Lon- 
don 

Carlsson, M., Hansteen, V, 2005, Chromospheric waves, in International Scientific Con- 
ference on Chromospheric and Coronal Magnetic Fields (ESA SP-596) on CD-Rom, 
(Eds.) D. Innes, A. Lagg, S. Solanki, D. Danesy 

Centeno, R., 2006, Investigacion de la propagacion de ondas en la atmosfera solar medi- 
ante espectropolarimetria en Helio I 10830 A, Ph.D. thesis, Universidad de La Laguna 

Centeno, R., Trujillo Bueno, J., Uitenbroek, H., CoUados, M., 2007, The influence of 
coronal EUV irradiance on the emission in the He I 10830 A and D3 multiplets, ArXiv 
e-prints, 712, ,Q712 . 22Q3, 

Chandrasekhar, S., 1960, Radiative Transfer, Dover Publications 

Choudhuri, A. R., Auffret, H., Priest, E. R., 1993, Implications of rapid footpoint motions 
of photospheric flux tubes for coronal heating, Sol. Phys., 143, 49-68 

CoUados, M., Lagg, A., Diaz Garcia, J. J., Hernandez Suarez, E., Lopez Lopez, R., Paez 
Maiia, E., Solanki, S. K., 2007, Tenerife Infrared Polarimeter II, in The Physics of 
Chromospheric Plasmas, (Eds.) P. Heinzel, I. Dorotovic, R. J. Rutten, vol. 368 of PASP 
Conf. Ser.,pp. 611-616 

de Boer, C. R., 1996, Noise filtering in solar speckle masking reconstructions, A&AS, 
120, 195-199 

De Moortel, I., Ireland, J., Hood, A. W., Walsh, R. W., 2002, The detection of 3& 5 min 
period oscillations in coronal loops, A&A, 387, 13-16 

De Pontieu, B., Erdelyi, R., James, S. P., 2004, Solar chromospheric spicules from the 
leakage of photospheric oscillations and flows. Nature, 430, 536-539 



96 



Bibliography 



De Pontieu, B., Mcintosh, S. W., Hansteen, V. H., Carlsson, M., Schrijver, C. J., Tarbell, 
T. D., Title, A. M., Shine, R. A., Suematsu, Y., Tsuneta, S., Katsukawa, Y., Ichimoto, 
K., Shimizu, T, Nagata, S., 2007, A Tale Of Two Spicules: The Impact of Spicules on 
the Magnetic Chromosphere, ArXiv e-prints, 710, Q71Q . 2934| 

Defouw, R. J., 1976, Wave propagation along a magnetic tube, ApJ, 209, 266-269 

Ferraro, V. C. A., Plumpton, A., 1966, An Introduction to Magneto-Fluid Mechanics, 
Oxford University Press, Oxford 

Fontenla, J. M., Avrett, E. H., Loeser, R., 1991, Energy balance in the solar transition 
region. 11 - effects of pressure and energy input on hydrostatic models, ApJ, 377, 712- 
725 

Fried, D. L., 1965, Statistics of a Geometric Representation of Wavefront Distortion, J. 
Opt. Soc. Am. (1917-1983), 55, 1427-1435 

Georgakilas, A. A., Alissandrakis, C. E., Zachariadis, T. G., 1990, Mass motions associ- 
ated with H-alpha active region arch structures, Sol. Phys., 129, 277-293 

Giovanelli, R., 1975, Wave systems in the chromosphere, Sol. Phys., 44, 299-314 

Grossmann-Doerth, U., Schmidt, W., 1992, Chromospheric fine structure revisited, A&A, 
264, 236-242 

Hansteen, V. H., De Pontieu, B., Rouppe van der Voort, L., van Noort, M., Carlsson, M., 
2006, Dynamic fibrils are driven by magnetoacoustic shocks, ApJ, 647, 73-76 

Heinzel, P., Dorotovic, I., Rutten, R. J. (Eds.), 2007, The Physics of Chromospheric Plas- 
mas, 368, ASP Conference Series 

Hildner, E., 1974, The formation of solar quiescent prominences by condensation, 
Sol. Phys., 35, 123-136 

Innes, D. E., Lagg, A., Solanki, S. A., Danesy, D. (Eds.), 2005, International Scientific 
Conference on Chromospheric and Coronal Magnetic Fields, vol. 596 of ESA Special 
Publication (CD-ROM) 

Janssen, K., 2003, Structure and dynamics of small scale magnetic fields in the solar 
atmosphere. Results of high resolution polarimetry and image reconstruction, Ph.D. 
thesis, Universitats-Sternwarte Gottingen 

Keller, C. U., von der Liihe, O., 1992, Solar speckle polarimetry, A&A, 261, 321-328 

Kippenhahn, R., MoUenhoff, C, 1975, Elementare Plasmaphysik, Bibliograpisches Insti- 
tut, Zurich 

Korff, D., 1973, Analysis of a Method for Obtaining Near Diffraction Limited Information 
in the Presence of Atmospheric Turbulence, J. Opt. Soc. Am. (1917-1983), 63, 971-980 

Koschinsky, M., Kneer, F, Hirzberger, J., 2001, Speckle spectro-polarimetry of solar 
magnetic structures, A&A, 365, 588-597 



97 



Bibliography 



Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., Tachikawa, S., Hashimoto, T., 
Minesugi, K., Ohnishi, A., Yamada, T., Tsuneta, S., Hara, H., Ichimoto, K., Suematsu, 
Y, Shimojo, M., Watanabe, T., Shimada, S., Davis, J. M., Hill, L. D., Owens, J. K., 
Title, A. M., Culhane, J. L., Harra, L. K., Doschek, G. A., Golub, L., 2007, The Hinode 
(Solar-B) Mission: An Overview, Sol. Phys., 243, 3-17 

Kukhianidze, V., Zaqarashvili, T. V., Khutsishvili, E., 2006, Observation of kink waves in 
solar spicules, A&A, 449, L35-L38, ar Xiv : astro- ph/Q6 Q2534| 

Labeyrie, A., 1970, Attainment of Diffraction Limited Resolution in Large Telescopes by 
Fourier Analysing Speckle Patterns in Star Images, A&A, 6, 85-87 

Lagg, A., 2007, Recent advances in measuring chromospheric magnetic fields in the He I 
10830 A line. Advances in Space Research, 39, 1734-1740 

Lagg, A., Woch, J., Krupp, N., Solanki, S. K., 2004, Retrieval of the full magnetic vector 
with the He I multiplet at 1083 nm. Maps of an emerging flux region, A&A, 414, 
1109-1120 

Lofdahl, M. G., 2002, Image Reconstruction from Incomplete Data II, in SPIE, (Eds.) 
P J. Bones, M. A. Fiddy, R. P Millane, vol. 4792, pp. 146-155 

Lofdahl, M. G., van Noort, M. J., Denker, C, 2007, Solar image restoration, in Mod- 
em solar facilities - advanced solar science, Gottingen September 27-29, 2006, (Eds.) 
F. Kneer, K. G. Puschmann, A. D. Wittmann, p. 1 19, Universitatsverlag Gottingen 

Martinez Pillet, V., CoUados, M., Sanchez Almeida, J., Gonzalez, V., Cruz-Lopez, A., 
Manescau, A., Joven, E., Paez, E., Diaz, J., Feeney, O., Sanchez, V., Scharmer, G., 
Soltau, D., 1999, LPSP & TIP: FuU Stokes Polarimeters for the Canary Islands Ob- 
servatories, in High Resolution Solar Physics: Theory, Observations, and Techniques, 
(Eds.) T. R. Rimmele, K. S. Balasubramaniam, R. R. Radick, vol. 183 of PASP Conf. 
Ser., pp. 264-272 

Mitalas, R., Sills, K. R., 1992, On the photon diffusion time scale for the sun, ApJ, 401, 
759-760 

Musielak, Z. E., Routh, S., Hammer, R., 2007, Cutoff-free Propagation of Torsional 
Alfven Waves along Thin Magnetic Flux Tubes, ApJ, 659, 650-654 

Neckel, H., 1999, Announcement spectral atlas of solar absolute disk-averaged and disk- 
center intensity from 3290 A to 125 10 A (Brault and Neckel, 1987) now available from 
Hamburg observatory anonymous FTP site, Sol. Phys., 184, 421-422 

Noll, R. J., 1976, Zernike polynomials and atmospheric turbulence, J. Opt. Soc. Am. 
(1917-1983), 66, 207-211 

Paxman, R. G., Seldin, J. H., Lofdahl, M. G., Scharmer, G. B., Keller, C. U., 1996, Evalu- 
ation of Phase-Diversity Techniques for Solar-Image Restoration, ApJ, 466, 1087-1099 

Priest, E. R., 1984, Solar Magnetohydrodynamics, Reidel, Dortrecht 



98 



Bibliography 



Puschmann, K. G., Sailer, M., 2006, Speckle reconstruction of photometric data observed 
with adaptive optics, A&A, 454, 1011-1019 

Puschmann, K. G., Kneer, R, Seelemann, T, Wittmann, A. D., 2006, The new Gottin- 
gen Fabry-Perot spectrometer for two-dimensional observations of the Sun, A&A, 451, 
1151-1158 

Robbrecht, E., Verwichte, E., Berghmans, D., Hochedez, J. R, Poedts, S., Nakariakov, 
V. M., 2001, Slow magnetoacoustic waves in coronal loops: EIT and TRACE, A&A, 
370, 591-601 

Roddier, N. A., 1990, Atmospheric wavefront simulation and Zemike polynomials, in 
Amplitude and Intensity Spatial Interferometry, (Ed.) J. B. Breckinridge, vol. 1237 of 
SPIE, pp. 668-679 

Rouppe van der Voort, L. H. M., De Pontieu, B., Hansteen, V. H., Carlsson, M., van 
Noort, M., 2007, Magnetoacoustic shocks as a driver of quiet-sun mottles, ApJ, 660, 
169-172 

Saha, S. K., 2002, Modern optical astronomy: technology and impact of interferometry. 
Reviews of Modern Physics, 74, 551-600, a rXiv : astro^ h/Q202 115] 

Sanchez- Andrade Nuno, B., Puschmann, K. G., Sanchez Cuberes, M., Blanco Rodriguez, 
J., Kneer, R, 2005, Chromospheric Dynamics of a Solar Active Region, in The Dy- 
namic Sun: Challenges for Theory and Observations, vol. 600 of ESA Special Publica- 
tion 

Sanchez- Andrade Nuno, B., Bello Gonzalez, N., Blanco Rodriguez, J., Kneer, R, 
Puschmann, K. G., 2007, Rast events and waves in solar Ha structures at high reso- 
lution, A&A, submitted 

Sanchez- Andrade Nuno, B., Centeno, R., Puschmann, K. G., Trujillo Bueno, J., Blanco 
Rodriguez, J., Kneer, R, 2007a, Spicule emission profiles observed in He I 10830 A, 
A&A, 472, 51-54, arXiv:Q7Q7.4421 

Sanchez- Andrade Nuno, B., Puschmann, K. G., Kneer, R, 2007b, Observations of a flar- 
ing active region in Hor, in Modern solar facilities - advanced solar science. Proceedings 
of a Workshop held at Gottingen, (Eds.) R. Kneer, K. G. Puschmann, A. D. Wittmann, 
pp. 273-276 

Secchi, P., 1877, Le Soleil, 2 

Soltau, D., 1985, The German 60-cm Vacuum Tower Telescope and Its Post-Rocus Ra- 
cilities, in Solar Physics and Interplanetary Travelling Phenomena, (Eds.) C. de Jager, 
B. Chen, pp. 1191-1198 

Spruit, H. C, 1981, Motion of magnetic flux tubes in the solar convection zone and chro- 
mosphere, A&A, 98, 155-160 

Spruit, H. C, 1982, Propagation speeds and acoustic damping of waves in magnetic flux 
tubes, Sol. Phys., 75, 3-17 



99 



Bibliography 



Sterling, A. C, 2000, Solar Spicules: A Review of Recent Models and Targets for Future 
Observations - (Invited Review), Sol. Phys., 196, 79-111 

Sterling, A. C, HoUweg, J. V., 1989, A rebound shock mechanism for solar fibrils, ApJ, 
343, 985-993 

Stix, M., 2002, The Sun. An Introduction, Springer, Berlin. Second Edition 

Tandberg-Hanssen, E., 1977, Prominences, in Illustrated Glossary for Solar and Solar- 
Terrestrial Physics, (Eds.) A. Bruzek, C. J. Durrant, vol. 69 of Astrophysics and Space 
Science Library, pp. 97-1 1 1 

Tomczyk, S., Mcintosh, S. W., Keil, S. L., Judge, P G., Schad, T, Seeley, D. H., Ed- 
mondson, J., 2007, Alfven Waves in the Solar Corona, Science, 317, 1192-1196 

Tothova, D., Innes, D. E., Solanki, S. K., 2007, Wavelet-based method for coronal loop 
oscillations analysis, in Modem solar facilities - advanced solar science. Proceedings 
of a Workshop held at Gottingen, (Eds.) F. Kneer, K. G. Puschmann, A. D. Wittmann, 
pp. 265-268 

Tripathi, D., Solanki, S. K., Mason, H. E., Webb, D. R, 2007, A bright coronal downflow 
seen in multi-wavelength observations: evidence of a bifurcating flux-rope?, A&A, 
472, 633-642 

Trujillo Bueno, J., Merenda, L., Centeno, R., CoUados, M., Landi Degl'Innocenti, E., 
2005, The Hanle and Zeeman Effects in Solar Spicules: A Novel Diagnostic Window 
on Chromospheric Magnetism, ApJ, 619, 191-194 

Tsiropoula, G., 2000, Physical parameters and flows along chromospheric penumbral 
fibrils, A&A, 357, 735-742 

Tsiropoula, G., Schmieder, B., 1997, Determination of physical parameters in dark mot- 
tles., A&A, 324, 1183-1189 

Tsiropoula, G., Tziotziou, K., 2004, The role of chromospheric mottles in the mass bal- 
ance and heating of the solar atmosphere, A&A, 424, 279-288 

Tsiropoula, G., Alissandrakis, C. E., Schmieder, B., 1993, The fine structure of a chro- 
mospheric rosette, A&A, 271, 574-586 

Tziotziou, K., 2007, Chromospheric Cloud-Model Inversion Techniques, in The Physics 
of Chromospheric Plasmas, (Eds.) P. Heinzel, I. Dorotovic, R. J. Rutten, vol. 368 of 
PASP Conf. Ser., pp. 217-237 

Tziotziou, K., Tsiropoula, G., Mein, P., 2003, On the nature of the chromospheric fine 
structure. I. Dynamics of dark mottles and grains, A&A, 402, 361-372 

Tziotziou, K., Tsiropoula, G., Mein, P., 2004, On the nature of the chromospheric fine 
structure. 11. Intensity and velocity oscillations of dark mottles and grains, A&A, 423, 
1133-1146 



100 



Bibliography 



van Noort, M., Rouppe van der Voort, L., Lofdahl, M. G., 2005, Solar Image Restora- 
tion By Use Of Multi-frame Blind De-convolution With Multiple Objects And Phase 
Diversity, Sol. Phys., 228, 191-215 

van Noort, M. J., Rouppe van der Voort, L. H. M., 2006, High-Resolution Observations 
of Fast Events in the Solar Chromosphere, ApJ, 648, 67-70 

Volkmer, R., Kneer, R, Bendlin, C, 1995, Short-period waves in small-scale magnetic 
flux tubes on the Sun., A&A, 304, 1^ 

von der Liihe, O., 1984, Estimating Fried's parameter from a time series of an arbitrary 
resolved object imaged through atmospheric turbulence, J. Opt. Soc. Am., 1, 510-519 

von der Liihe, O., Sohau, D., Berkefeld, T., Schelenz, T, 2003, KAOS: Adaptive optics 
system for the Vacuum Tower Telescope at Teide Observatory, in Innovative Telescopes 
and Instrumentation for Solar Astrophysics., (Eds.) S. L. Keil, S. V. Avakyan, vol. 4853 
of SPIE Conf., pp. 187-193 

von Uexkiill, M., Kneer, R, Mattig, W., 1983, The chromosphere above sunspot umbrae, 
rv - Frequency analysis of umbral oscillations, A&A, 123, 263-270 

Weigelt, G., Wimitzer, B., 1983, Image reconstruction by the speckle-masking method. 
Optics Letters, 8, 389-391 

Weigelt, G. P., 1977, Modified astronomical speckle interferometry 'speckle masking'. 
Optics Communications, 21, 55-59 

Wentzel, D. G., 1979, Hydromagnetic surface waves on cylindrical fluxtubes, A&A, 76, 
20-23 

Wikipedia, t. f. e.. Stokes parameters, (http://en.wikipedia.org/wiki/Stokes_parameters) 
Wikipedia, t. f. e.. Sun, (http://en.wikipedia.org/wiki/sun) 

Wilhelm, K., 2000, Solar spicules and macrospicules observed by SUMER, A&A, 360, 
351-362 

Wittmann, A., 1969, Some properties of umbral flashes, Sol. Phys., 7, 366-369 

Yi, Z., Darvann, T., Molowny-Horas, R., 1992, Software for Solar Image Processing - 
Proceedings from lest Mini Workshop, Tech. Rep. 56, Lest foundation. Institute of 
Theoretical Astrophysics, University of Oslo 



101 



Publications 



Refereed papers 

1. B. Sanchez- Andrade Nuno, R. Centeno, K. G. Puschmann, J. Trujillo Bueno, 

o 

J. Blanco Rodriguez, and F. Kneer. Spicule emission profiles observed in He 1 10830 A. 
A&A, 472:L51-L54, Sept. 2007. 

2. B. Sanchez-Andrade Nuno, N. Bello Gonzalez, J. Blanco Rodriguez, F. Kneer, 
and K. G. Puschmann. Fast events and waves in an active region of the Sun ob- 
served in Ha with high spatial resolution. A&A, submitted, Dec. 2007. 

3. J. Blanco Rodriguez, O. V. Okunev, K. G. Puschmann, F. Kneer, and B. Sanchez- 
Andrade Nuno. On the properties of faculae at the poles of the Sun. A&A, 
474:251-259, Oct. 2007. 

Conference contributions 

4. B. Sanchez-Andrade Nunc, R. Centeno, K. G. Puschmann, J. Trujillo Bueno, and 
F. Kneer. OfF-limb spectroscopy of the He I 10830 A multiplet: observations vs. 
modelling. In F. Kneer, K. G. Puschmann, and A. D. Wittmann, editors. Modem 
solar facilities - advanced solar science, pages 177-180, 2007. 

5. B. Sanchez-Andrade Nuno, K. G. Puschmann, and F. Kneer. Observations of a 
flaring active region in Hor. In F. Kneer, K. G. Puschmann, and A. D. Wittmann, 
editors. Modem solar facilities - advanced solar science, pages 273-276, 2007. 

6. B. Sanchez-Andrade Nuno, K. G. Puschmann, M. Sanchez Cuberes, J. Blanco 
Rodriguez, and F. Kneer. Analysis of a Wide Chromospheric Active Region. In 
D. E. Innes, A. Lagg, S. A. Solanki, and D. Danesy, editors, Chromospheric and 
Coronal Magnetic Fields, volume 596 of ESA Special Publication, Nov. 2005. 

7. B. Sanchez-Andrade Nuno, K. G. Puschmann, M. Sanchez Cuberes, J. Blanco 
Rodriguez, and F. Kneer. Chromospheric Dynamics of a Solar Active Region. In 
The Dynamic Sun: Challenges for Theory and Observations, volume 600 of ESA 
Special Publication, Dec. 2005. 

8. B. Sanchez-Andrade Nuno. Study case: Solar Science Communication. In 
L. Lindberg Christensen, and M. Zoulias, editors. Communicating Astronomy with 
the Public 2007. An lAU/Nat. Obs. of Athens/Eugenides Foundation Conference, 
Oct. 2007. 



103 



9. J. Blanco Rodriguez, B. Sanchez-Andrade Nuno, K. G. Puschmann, and F. Kneer. 
Study of Polar Faculae. In The Dynamic Sun: Challenges for Theory and Obser- 
vations, volume 600 of ESA Special Publication, Dec. 2005. 

10. J. Blanco Rodriguez, B. Sanchez-Andrade Nuno, K. G. Puschmann, and F. Kneer. 
Study of Polar Faculae. In D. E. Innes, A. Lagg, S. A. Solanki, and D. Danesy, 
editors, Chromospheric and Coronal Magnetic Fields, volume 596 of ESA Special 
Publication, Nov. 2005. 

11. F. Kneer, K. G. Puschmann, J. Blanco Rodriguez, B. Sanchez-Andrade Nuno, 

and A. D. Wittmann. Magnetic Structures on the Sun: Observations with the 
New "GOTTINGEN" Two-Dimensional Spectrometer on Tenerife. In D. E. Innes, 
A. Lagg, and S. A. Solanki, editors, Chromospheric and Coronal Magnetic Fields, 
volume 596 of ESA Special Publication, Nov. 2005. 

12. L. Valdivielso Casas, N. Belle Gonzalez, K. G. Puschmann, B. Sanchez-Andrade 
Nuno, and F. Kneer. Analysis of Polarimetric Sunspot Data from Tesos/vtt/tenerife. 
In D. E. Innes, A. Lagg, S. A. Solanki, and D. Danesy, editors, Chromospheric and 
Coronal Magnetic Fields, volume 596 of ESA Special Publication, Nov. 2005. 

13. C. Denker, A.P. Verdoni, F. Woger, A. Tritschler, T.R. Rimmele, F. Kneer, K. Rear- 
don, B. Sanchez-Andrade Nuno, I. Dominguez Cerdena and K.G. Puschmann 

Speckle Interferometry of Solar Adaptive Optics Imagery DFG-NSF Astrophysics 
Research Conference "Advanced Photonics in Application to Astrophysical Prob- 
lems", June 2007 



Acknowledgements 



I would like to acknowledge here the contributions from many people to the successful 
completion of this work. I want to express my gratitude to all those who helped me, from 
my supervisor this three last years to my teachers at my first university who encouraged 
me to start this adventure in astrophysics. During all this time I have been also surrounded 
and supported by many people to whom I want to thank within these lines. They all made 
this experience of doing a PhD, professionally and personally, one of the best times of my 
Ufe. 

My supervisor of this research work was Prof. Dr. Franz Kneer, or better /ranz. He 
was the first person I saw when I came to the train station back three years ago (along with 
Julian). He was the guide to all my work, the observations, the interpretation, . . . basically, 
my formation as solar physicist. The amount of knowledge I learnt from him cannot be 
measured by any means, and for that I am professionally extremely grateful. In a personal 
way he was also very kind and helped me always with good advices whenever I needed 
it. He understood me perfectly when I had problems and also encouraged me to travel as 
much as I wanted (which was not a little). Working with him this time made the whole 
experience of PhD leaps better. Ich kann nicht diese Arbeit ohne Ihn zu vorstellen. Danke 
Franz. 

Above all, I am and will be always grateful to my reference in life, my parents Con- 
chita Nuno Lopez and Julio Sanchez- Andrade Fernandez, and my sister Deva Sanchez 
Nuno. They always supported me, although a bit far in physical distance during this last 
years. When I was 5 years old I was going home with my mother from the playground 
when I saw a poster announcing a conference about starts. I asked my mother to read it 
loud for me, and to her surprise I had to explained her: "Don't you yet know that I want 
to be a researcher ofstarsT. The day after she gave me a children book about stars and 
explained me that if I wanted to be so I should know that they are called astronomers. 
And that's how it all began. 

My work was supported by two institutions: the Max Planck Institute fur Sonnen- 
systemforschung (MPS) granted me the fellowship and the Institut fur Astrophysik Got- 
tingen (TAG) provided me the facilities to work with Franz at Gottingen and support 
for the observations at Tenerife. Also, being part of the International Max Planck Re- 
search School on Physical Processes in the Solar System and Beyond at the Universities 
of Braunschweig and Gottingen (IMPRS) I could attend seminars, retreats and courses on 
various astrophysical topics. I am very grateful to these institutions for providing such a 
broad curricula. I would like also to thank the coordinator of the MPS and IMPRS, Dr. 
Dieter Schmitt. 

The Vacuum Tower Telescope used for the observations is operated by the Kiepenheuer- 
Institut fiir Sonnenphysik, Freiburg, at the Spanish Observatorio del Teide of the Instituto 



105 



Acknowledgements 



de Astrofisica de Canarias. 

Here at the solar group of TAG we had a really stimulating environment with long 
discussions: Franz, Klaus, Markus, Nazaret, Julian and the various Erasmus people that 
went by (like Luisa, Manu, . . . )• We worked in the beginning at the beautiful historical 
Sternwarte and then at the Faculty of Physics. 1 would like to thank specially Klaus for 
all the help and scientific discussions during all the time he was here. 

Many other professional colleagues contributed to this work. Their input was ex- 
tremely helpful. A short list with few names would include: (from the MPS) Andreas 
Lagg, Vasili Zakarov, (from the lAG) the solar group. Axel Wittman, Volker Bothmer, 
(from the lAC) Basilio Ruiz Cobos, Manolo CoUados, Javier Trujillo, Rebecca Centeno. 
Thanks to Nurol Al and H. Schleicher for their computing codes. For the Blind Decon- 
volution section I got much help from Jaime and Mats Lofdahl at the Institute for Solar 
Physics, Stockholm. 

From Asturias to here there is long way, specially passing through Canary Island. 
Without the advices and guidance from many people I would have been lost. In these 
sense I specially thank Cristina, Juanjo, Basilio and Fernando. 

Throughout all these years I had the luck to be surrounded by the best companion of 
friends. I would like to mention lastly some few friends that helped me to keep my feet on 
the ground. Here in Gottingen Thorsten, Gonzaga, Klaus, Miguel, Cristiano, Vladi, Cathi, 
Niko, Lome, Diego, Iria, Teresa, Nora, Carlos, Olga, Julian, Nazaret, Benoit, Markus, 
Manu and Luisa. Thorsten in particular was not only a good friend but also my flat-mate, 
climbing-mate, snowboard-mate and of course party-mate. My new flat-mate, Richard, 
had to stand me this last months with my thesis mood, though. Thanks a lot!. Mark, 
Lucas, Clementina, Pedro, Emre or Martin, all the guys from Lindau, always warmly 
welcomed his lost member of the family emigrated to the civilization. 

During my short intense period in Tenerife I made good friends: Onti, Borja, Bendinat, 
Adriana, Manu, Miguel. Some time later I met Rebecca, with whom I also worked these 
years. 

Whenever a crazy trip came up I had AEGEE people coming along. They would join 
me to any place in Europe. With them I had some of my best times: Luis, David, Neila, 
Saioa, Javiero, Marti, Marta and recently of course Carol. All those friends with crazy 
names: Annia, Annamary, Bir, Konstantina, Marios or Lutzn. I am sure I'll meet you 
again, somewhere. 

Since I left Asturias, I have always missed the green and astounding landscapes, but 
also my old friends: Roberto, Raul, Estela, Flaci, Nacho, Cova, Lisa, David and of course 
my big and warm family. 

Por ultimo me gustaria terminar esta seccion, y con ello el fin de este trabajo, con la 
memoria de tres personas a quienes siempre echare de menos. En mi primer ano en Ale- 
mania tuve que despedir a mi abuela matema y el segundo a mi abuela paterna. Aunque 
sea ley de vida no deja de ser doloroso. Este ultimo ano, a mi primo Abraham, a quien 
tuve la suerte de ver una ultima vez. Con ellos he aprendido la leccion mas importante: 

Disfruten de la vida. 

Gottingen, January 2008 



106 



Lebenslauf 



Name: 

Geburtsdatum: 
Geburstort: 
Familienstand: 
Eltem: 



Bruno Sanchez-Andrade Nuno 
6. Mai 1981 

Oviedo, Asturias, Spanien 
Ledig 

Julio Miguel Sanchez-Andrade Fernandez 
Concepcion Nuno Lopez 



Staatsangehorigkeit: Spanisch 



Schulbildung: 



Studium: 



Promotion: 



September 1987 - Juni 1995 Grundschule am Colegio 

Publico "Rio Piles" in Gijon, 
Asturias 

September 1995 - Juni 1999 Weiter fiihrende Schule am /. 

B. "Rio Piles" in Gijon, As- 
turias 

September 1999 - September 2003 Physikalische Fakultat der Uni- 

versitat Oviedo, Asturias 
September 2003 - September 2004 Physikalische Fakultat der Uni- 

versitat La Laguna, Teneriffa 
(Fachrichtung Astrophysik) 



Januar 2005 - Januar 2008 



Januar 2005 - Januar 2008 



Promotion an der Universitats- 
Sternwarte Gottingen (seit Juni 
2005 Institut fiir Astrophysik 
Gottingen, lAG) 
Stipendium des Max-Planck- 
Instituts fiir Sonnensystem- 
forschung 



107 



